arXiv:1507.06937v2 [hep-ph] 6 Dec 2015 


ICAS 02/15 
TIF-UNIMI-2015-10 
ZU-TH 24/15 


Vector boson production at hadron colliders: 
transverse-momentum resummation 
and leptonic decay 


Stefano Catani(“\ Daniel de Florian*^^\ 

Giancarlo Ferrera*^^) and Massimiliano Grazzini*^'^0 

(“)INFN, Sezione di Firenze and Dipartimento di Fisica e Astronomia, 
Universita di Firenze, 1-50019 Sesto Fiorentino, Florence, Italy 
*'^^Departamento de Fisica, FCEYN, Universidad de Bnenos Aires, 

(1428) Pabellon 1 Cindad Univer sit aria. Capital Federal, Argentina 
and International Center for Advanced Stndies (ICAS), 

UNSAM, Campus Miguelete, 25 de Mayo y Francia, 

1650 Buenos Aires, Argentina 
Dipartimento di Fisica, Universita di Milano and 
INFN, Sezione di Milano, 1-20133 Milan, Italy 
{'i)piiysik-lnstitut, Universitat Zurich, CH-8057 Zurich, Switzerland 

Abstract 

We consider the transverse-momentum [qt) distribution of Drell-Yan lepton pairs 
produced, via W and Zj'y* decay, in hadronic collisions. At small values of we 
resum the logarithmically-enhanced perturbative QCD contributions up to next-to- 
next-to-leading logarithmic accuracy. Resummed results are consistently combined 
with the known (P(a|) fixed-order results at intermediate and large values of qr- Our 
calculation includes the leptonic decay of the vector boson with the corresponding 
spin correlations, the hnite-width effects and the full dependence on the final-state 
lepton(s) kinematics. The computation is encoded in the numerical program DYRes, 
which allows the user to apply arbitrary kinematical cuts on the final-state leptons and 
to compute the corresponding distributions in the form of bin histograms. We present 
a comparison of our results with some of the available LHC data. The inclusion of the 
leptonic decay in the resummed calculation requires a theoretical discussion on the g-r 
recoil due to the transverse momentum of the produced vector boson. We present a 
g-r recoil procedure that is directly applicable to qt resummed calculations for generic 
production processes of high-mass systems in hadron collisions. 
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1 Introduction 


The production of high-mass lepton pairs through the Drell-Yan (DY) mechanism [T| is a bench¬ 
mark hard-scattering process at hadron colliders. It provides important tests of the Standard 
Model (SM) with both precise measurements of its fundamental parameters and, at the same 
time, stringent constraints on new physics. 

It is thus a major task to achieve accurate theoretical predictions for the DY production cross 
section and related kinematical distributions. This requires, in particular, the evaluation of QCD 
radiative corrections, which can be perturbatively computed as power series expansion in the 
strong coupling as- The total cross section [2] and the rapidity distribution |3] of the vector 
boson are known up to the next-to-next-to-leading order (NNLO) in perturbative QCD. Two 
independent fully exclusive NNLO calculations, which include the leptonic decay of the vector 
boson, have been performed ilElE]. Electroweak (EW) radiative corrections are also available 
for both W [7] and Zj^* [8] production. Mixed QCD-EW corrections have been considered in 
Refs. |9l[ini[IIl[l2]. 

A particularly relevant observable is the transverse-momentum (g^) distribution of the vector 
boson. To obtain a precise measurement of the W mass it is important to have accurate theoretical 
calculations of the W and Z bosons qt spectra. In the large-g^ region ~ my), where the 
transverse momentum is of the order of the vector boson mass my, QCD corrections are known 
up to 0(ag) [131HH US], and these results were extended in Refs. [Ml HT] with the inclusion of the 
dependence on the leptonic decay variables. Very recently the fully exclusive 0{a^) computation 
of vector boson production in association with a jet has been performed in Ref. [18] (in the case 
of W production) and Ref. [T9| (in the case of Zj'y* production). 

The bulk of the vector boson cross section is produced in the small-gr region (g^ <C my), 
where the reliability of the fixed-order expansion is spoiled by the presence of large logarithmic 
corrections, ag {rriy/q^) \n^{my/q^) (with 0 < m < 2n —1), of soft and collinear origin. To obtain 
reliable predictions, these logarithmically-enhanced terms have to be evaluated and systematically 
resummed to all orders in perturbation theory [20] [29]. In recent years, the resummation of small- 
g-r logarithms has been reformulated [30]-[38] by using Soft Collinear Effective Theory (SCET) 
methods and transverse-momentum dependent (TMD) factorization. 

The resummed and fixed-order calculations, which are valid at small and large values of qr, 
respectively, can be consistently matched at intermediate values of qx to achieve a uniform theo¬ 
retical accuracy for the entire range of transverse momenta. 

In this paper we compute the vector boson transverse-momentum distribution [SHI HD] by 
using the resummation formalism proposed in Refs. [23123 EZ], which can be applied to a generic 
process in which a high-mass system of non strongly-interacting particles is produced in hadronic 
collisions [1I11261127], [l2]-[52]. Other phenomenological studies of the DY qx distribution, which 
combine resummed and hxed-order perturbative results at different levels of theoretical accuracy, 
can be found in Refs. [S3]^[HH]- Within the studies in Refs. [M]^[n2]) the kinematical dependence 
on the momenta of the final-state leptons is considered only in the RESBOS calculation [5511551155] 
and in the calculations of Refs. EH and [57] . 

Hadron collider experiments can directly measure only the decay products of vector bosons in 
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finite kinematical regions. Therefore, it is important to include the vector boson leptonic decay in 
the theoretical calculations, by retaining the kinematics of the hnal-state leptons. In this way it is 
possible to obtain predictions for the transverse-momentum distribution of the measured leptons. 
This is specially relevant in the case of W production where, because of the hnal-state neutrino, 
the transverse momentum of the vector boson can only be reconstructed through a measure of 
the hadronic recoil. Moreover, in both cases of W and Z production, the inclusion of the leptonic 
decay allows one to apply kinematical selection cuts, thus providing a more realistic simulation of 
the actual experimental analysis. 

In Ref. [39l ITO] we have presented a resummed computation of the transverse-momentum 
spectrum for Z/y* production at Tevatron energies. We have combined resummation at the next- 
to-next-to-leading logarithmic (NNLL) accuracy in the small-gr region with the hxed-order results 
at 0{a‘g) in the large-g^ region. This leads to a calculation with uniform theoretical accuracy 
from small to intermediate values of qt- In particular, the integral over the range 0 < gr < gTmax 
(gumax is a generic upper limit in the small-gr region) of the qr distribution includes the complete 
perturbative terms up to NNLO. Moreover, at large values of qt the calculation implements a uni- 
tarity constraint that guarantees to exactly reproduce the NNLO value of the total cross section 
after integration over qt- In this paper we extend the NNLL-I-NNLO calculation of Ref. m to 
W boson production, and we include the leptonic decay of the vector boson with the correspond¬ 
ing spin correlations. The spin of the vector boson dynamically correlates the decaying lepton 
momenta with the transverse momentum acquired by the vector boson through its production 
mechanism. Therefore, the inclusion of the full dependence on the lepton decay variables in the 
resummed calculation requires a theoretical discussion on the treatment of the qt recoil due to 
the transverse momentum of the vector boson. We treat the qt recoil by introducing a general 
procedure that is directly applicable to qt resummed calculations for generic production processes 
of high-mass systems in hadron collisions. The calculation presented in this paper parallels the one 
performed in Ref. [13] for the case of SM Higgs boson production, with the non-trivial additional 
complication of dealing with the spin correlations that are absent in the Higgs boson case. Our 
vector boson computation is implemented in the numerical code DYRes, which allows the user 
to apply arbitrary kinematical cuts on the hnal-state leptons and to compute the corresponding 
relevant distributions in form of bin histograms. The code DYRes is publicly available and it can 
be downloaded from the URL address http: //pcteserver .mi . infn. it/~ferrera/dyres .html . 

The paper is organized as follows. In Sect. [2] we briehy review the resummation formalism of 
Refs. [23, lini |27j , and we discuss the main features of qt resummation for the DY process with full 
dependence on hnal-state lepton variables. In Sect. [3] we present our quantitative results for vector 
boson production at LHC energies. Section 13.11 is devoted to the qt spectrum of the vector boson 
after integration over the hnal-state leptons. We present results at diherent orders of logarithmic 
accuracy, we study the corresponding dependence on scale variations, and we briehy comment on 
uncertainties due to parton densities and on non-perturbative ehects. In Sect. 13.21 we compare 
our numerical results for Z/y* and W production with some of the available LHC data, and we 
also study the impact of transverse-momentum resummation on lepton kinematical variables. In 
Sect. Hlwe summarize our results. The Appendix presents a detailed discussion of qx recoil and of 
its implementation. 
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2 Transverse-momentum resummation 


In this Section we briefly recall the main features of the transverse-momentum resummation 
formalism that we use in this paper. A more detailed discussion of the resummation formalism 
can be found in Refs. [25| l26| \T7\ [29]. In Ref. [10] we have considered NNLL resummation for 
the qt distribution of the vector boson after integration over the kinematical variables of the 
decaying leptons and the rapidity of the vector boson. In this paper we extend the results of 
Ref. [ID] to include the entire kinematical dependence on the final-state leptons. The presentation 
in this Section parallels that of Sect. 2 in Ref. [40] and, in particular, we highlight the main 
differences that arise in the treatment of the rapidity of the vector boson and, especially, of the 
lepton kinematics. 


We consider the inclusive hard-scattering process 


hi{Pi) + h 2 {P 2 ) ^ ^3(^3) + ^4(^4) + W, (1) 


where the collision of the two hadrons hi and ^2 with momenta Pf and P 2 produces the vector 
boson V {y = , W~ 1 Z and/or 7 *) with total momentum which subsequently decays in the 

lepton pair and X denotes the accompanying final-state radiation. We consider high values 
of the invariant mass M of the lepton pair (in general, M differs from the on-shell mass mv of the 
vector boson V), and we treat the colliding hadrons and the leptons in the massless approximation 
{Pf = P 2 = P 3 = P 4 = 0) throughout the paper. In a reference frame where the colliding hadrons 
are back-to-back, the momentum is fully specified by the invariant mass M (M^ = g^), the 
two-dimensional transverse-momentum vector (with magnitude qt = and azimuthal angle 
and the rapidity y {y = \ \n of the vector boson. Analogously the momentum of the 

lepton Ij (j = 3,4) is specified by the lepton rapidity pj and transverse momentum pxj- 


The kinematics of the lepton pair is completely specified by six independent variables (e.g., 
the three-momenta of the two leptons). For our purposes, it is convenient to use the vector boson 
momentum g^ to select four independent variables. Therefore, the final-state lepton kinematics 
is fully determined by the vector boson momentum g^ = P 3 + P 4 and by two additional and 
independent variables that specify the angular distribution of the leptons with respect to the 
vector boson momentum g^. We generically denote these two additional kinematical variables as 
Q, = These two independent variables can be chosen in different ways. For instance, 

we can use longitudinally boost invariant variables such as the rapidity difference Ps — y and the 
azimuthal angle 03 (or the azimuthal angle difference 03 — 0 qj,) of the lepton I 3 and the vector 
boson in the hadronic back-to-back reference frame. Alternatively, we can use the polar and 
azimuthal angles { 6 ', 0 '} of one lepton in a properly specified rest frame of the vector boson (such 
as, for instance, the Collins-Soper rest frame EDI). Independently of the actual specification of 
the variables 12 , the most general fully-differential hadronic cross section is expressed in terms of 
the sixfold differential distribution 


de^h\h2 — 

d? Qt dM"^ dp dCl 


{qT,M,y,s, 12 ) , 


( 2 ) 


where s = {Pi - 1 -^ 2 )^ = 2Pi ■ P 2 is the square of the hadronic centre-of-mass energy. Obviously, 
the differential distribution also depends on the EW parameters (including the mass my of the 
vector boson V)\ unless otherwise specified, this dependence is not explicitly denoted throughout 
the paper. 
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The differential hadronic cross section can be written as 


da, 


hlh2—>-l3l4 


d? Qt dM'^ dy dil 


[qT,M,y,s,ft) = V / dxi [ 
ai.a2>^0 Jo 


da. 


X 


aia2^/3/4 


Qt dM‘^ dy dfl 


dX2 fai/hi (^1 ) /^_f) fa2lh2 (^2 ) /^_f) 


[qT,M,y,s,n;as{nR),HR,HF), ( 3 ) 


where fa/h{x, (a = qj, qf, g) are the parton densities of the colliding hadron h at the factoriza¬ 
tion scale /ii?, ddaja 2 _j.z 34 are the differential partonic cross sections, s = X\X 2 S is the square of the 
partonic centre-of-mass energy, y = y — \n. \/x\lx 2 is the vector boson rapidity with respect to 
the colliding partons, and /ir is the renormalization scale. Note that the partonic cross sections do 
not have any explicit dependence on hadronic kinematical variables, since the leptonic variables 
n are specihed with respect to q^. The partonic cross section dda^a 2 ^hi 4 is computable in QCD 
perturbation theory as a power series expansion in the QCD coupling as- 


In the region where qT ~ M, the perturbative expansion of the partonic cross section starts at 
0{as)- In this region the value of the auxiliary scales yp and yu can be chosen to be of the order 
of M, and the QCD perturbative series is controlled by a small expansion parameter as{M‘^). 
Therefore, hxed-order calculations of the partonic cross section are theoretically justihed. The 
QCD radiative corrections are known analytically up to 0(a|) after integration over the lepton 
angular variables dans] and with the inclusion of the full dependence on these angular variables 
[min]. The numerical results at C?(q!|) can be obtained also from the fully-exclusive calculations 
of Refs, [a El E]. Results at 0{ag) can be derived from the recent numerical computations of 
W -I-jet production [TS] and -|-jet production [T^j. 

In the small qp region (qp-^M), the perturbative computation of the partonic cross section 
starts at 0{a^g) through the leading-order (LO) EW process qfqy ^ R of quark-antiquark an¬ 
nihilation. In this region, the QCD radiative corrections are known up to NNLO (i.e., (9(a|)) 
in analytic form m by neglecting corrections of 0{qp/M) (these corrections can directly be ex¬ 
tracted from Refs. [la la ESI nil). The complete (i.e., by including corrections of 0{qp/M)) 
NNLO result can be obtained from the numerical computations of Refs. BE]. However, in the 
small qp region the convergence of the hxed-order perturbative expansion is spoiled by the pres¬ 
ence of powers of large logarithmic terms, Og {M'^/qp) ln'"(M^/g|.) (with 0 < m < 2n — 1). In 
particular, these terms become singular in the limit qp 0. To obtain reliable predictions these 
terms have to be resummed to all orders. 


Within our formalism, the resummation is performed at the level of the partonic cross section, 
which is decomposed as follows: 



a.ia 2 — 



(res.) 


+ 



(fin.) 


( 4 ) 


Here we have introduced a shorthand notation: the symbol 


da 


aia 2 —>-/3U 


denotes the multidiffer¬ 


ential partonic cross section that appears as the last factor in the right-hand side of Eq. ([3]). The 
hrst term, on the right-hand side of Eq. (jl]) is the resummed component. It contains all 

the logarithmically-enhanced contributions (at small qp) that have to be resummed to all orders 
in as- The second term, the finite component dd^^^'\ is free of such contributions and thus it 
can be evaluated at fixed order in perturbation theory. Note that part of the non-singular (i.e.. 
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not logarithmically-enhanced) contributions can also be included in and we comment later 

about this point. 


The resummation of the logarithmic contributions has to be carried out in the impact pa¬ 
rameter (b) space [20l [21], [22l (53] [23] to fulfil the important constraint of transverse-momentum 
conservation for inclusive multiparton radiation. The impact parameter b is the conjugate variable 
to Qt through a Fourier transformation. The small-gr region [qt -C M) corresponds to the large-6 
region [bM ^ 1) and the logarithmic terms ln(M^/g|,) become large logarithmic contributions 
ln(M^6^) in 6 space. The resummed component of the cross section is then obtained by performing 
the inverse Fourier transformation (or the Bessel transformation in Eq. ([6|)) from b space to Qt 
space. The resummed component of the partonic cross section in Eq. (jTj) can be expressed as 


da. 


(res.) 

aia 2 — 




dd. 


bl,b2=<]f,qfi 


( 0 ) 

^ 1 ^ 2—>-^34 

dCl 


f^aia2,6ifc2-)-v(5'T; M, y, S] /i^, , 


( 5 ) 


where 

f^aia 2 ,bi 62 —^ Vi ^Si ^'Ri 



b Joibgr) yVa^a 2 Mb 2 ^v{b, M, y, s; as, yp, yl) , 


( 6 ) 


and Jo{x) is the Oth-order Bessel function. The factor in the right-hand side of Eq. ([5|) 

is the Born level differential cross section for the partonic subprocess qjqf —)■ E —)■ of quark- 
antiquark annihilation, where the quark flavours / and /' can be either different (if V = IF^) or 
equal (if V = Z,'y*). This factor is of purely EW origin, and it completely encodes the dependence 
on the lepton kinematical variables fl. We postpone more detailed comments on dd^^^ (see Eq. (na 
and the discussion therein). The QCD radiative corrections and their associated dependence on 
ln(M^6^) are embodied in the resummed factor yVa.^a 2 ,bib 2 ^v, which depends on the produced 
vector boson V but it is independent of the decay leptons (in particular, it does not depend on 
fi). The integrand W in Eq. ([6D depends on b^ = 6^ and the inverse Fourier transformation 
is recast in terms of the Bessel transformation through the integration over the azimuthal angle 
of b. Note that the resummation factor Waia 2 ,bib 2 ^v depends on q^ and it does not contain 
any dependence on the azimuthal angle (pq^, of qr- This azimuthal independence is a feature of 
transverse-momentum resummation [23| for the production processes of colourless systems (such 
as vector bosons) through quark-antiquark annihilation. In contrast, logarithmically-enhanced 
azimuthal correlations enter transverse-momentum resummation for processes initiated by gluon- 
gluon fusion [28] (such as Higgs boson production) and for production of systems that carry colour 
charges (such as heavy quarks) [7T] through either quark-antiquark annihilation or gluon-gluon 
fusion. 


The all-order resummation structure of 'Waia 2 Mb 2 ^v in Eq. (|2]) can be organized in exponential 
form [26] [27]. The exponentiated structure is directly evident by considering the ‘double’ (W, N 2 ) 
Mellin moments M) of the function Wvib, M,y, s) with respect to the variables zi = 

e+i'M/v^ and Z 2 = e~'^M/\/I at hxed M. We have0 

wl^^’^^\b,M;as{yl),yl,yl) = n^y^’^^^{M;as{yl),M/yR,M/yp,M/Q) 

X exp{g^^^’^^\as{yR),L]M/yR,M/Q)} , (7) 

tPor the sake of simplicity, in this presentation we omit the explicit dependence on the parton indices 
{0102,6162}. This simplified notation applies to the case of a sole parton species or, more precisely, to flavour 
non-singlet partonic channels (see Refs. [251 [ 27 ] for the general case). 
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where the dependence on b (and on the large logarithm ln(M^ 6 ^)) is denoted by dehning and 
introducing the logarithmic expansion parameter L = ln((5^&^/feo + 1) with bo = {'Je = 

0.5772... is the Euler number). The scale Q ~ M, named resummation scale [H], which appears 
in the right-hand side of Eq. ([7]), parametrizes the arbitrariness in the resummation procedure. 
Although does not depend on Q when evaluated to all perturbative orders, its explicit 

dependence on Q occurs when it is computed by truncation of the resummed expression at some 
level of logarithmic accuracy (see Eq. (|H])). Variations of Q around M can thus be used to estimate 
the size of yet uncalculated higher-order logarithmic contributions. 

The contribution exp{^^'^i’'^ 2 )| the right-hand side of Eq. ([7]) includes the Sudakov form 
factor and collinear-evolution terms. This contribution (which does not depend on the factorization 
scale /ip) is universal (i.e. process independent), namely, it is independent on the produced 
vector boson V and, more generally, it occurs in transverse-momentum resummation for all the 
processes that are initiated by quark-antiquark annihilation at the LO level. The generalized 
form factor exp{^*'^i’^ 2 ^} contains all the terms that order-by-order in as are logarithmically 
divergent as 6 —>■ cxo (or, equivalently, as qt —>■ 0). The all-order expression of the form factor 
can be systematically expanded in terms of functions g^^\asL) of the resummation parameter 
(each function g^^\asL) resums terms a^L^ and it is dehned such that (7*'^^(0) = 0 ). 
The resummed logarithmic expansion of in powers of 05 (p|.) reads 




= Lg^^\asinE)L) + M/M/Q) 

+ + , ( 8 ) 

TT 


where the term L g^^'^ collects the leading logarithmic (LL) contributions, the function g^‘^^ includes 
the next-to-leading logarithmic (NLL) contributions |23], g^^'> controls the NNLL terms [72l |32] 
and so forth. The function depends on the specific process of vector boson production and 

it is due to hard-virtual and collinear contributions. This function does not depend on the impact 
parameter b (it includes all the perturbative contributions to that behave as constants 

in the limit b 00 ) and, therefore, it can be expanded in powers of as = asifiji) as 

71 \ 71 / 

The next-to-leading order (NLO) term 7 {d)(Wi,Af 2 ) j^nown since a long time |73], and the NNLO 
term '^^AWi,Af 2 ) been obtained more recently by two independent calculations in Refs. [71] 
and m The explicit form of the functions -^d)(Wi,Af 2 ) and, in particular, their 

dependence on the Mellin moment indices (iVi, V 2 ) can be found in Ref. [26| and in Appendix A 
of Ref. [27| . 

Incidentally, we recall that the generalized form factor exp{^} is known up to NNLL accuracy 
also for processes initiated by the gluon fusion mechanism uni mi ESI Eg, and that the Cl(a|) 
collinear coefficients (which contribute to the NNLO term in Eq. dOj)) are also known for all possible 
partonic channels [THl EH EH ESI [29] . Owing to the universality structure of transverse-momentum 
resummation, these results and those for the qq annihilation channel (which contribute to vector 
boson production) can be directly implemented in resummed calculations for production processes 
of generic high-mass systems. 
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The finite component in Eq. (j3]) has to be evaluated starting from the usual fixed- 

order perturbative truncation of the partonic cross section and subtracting the expansion of the 
resummed part at the same perturbative order. We have 



(fin.) 

aia 2 — 


f.o. 



aia 2 — 


f.o. 



(res.) 

aia2—>-/3/4 


f.o. 


( 10 ) 


where the subscript f.o. denotes the perturbative truncation at the order f.o. (NLO, NNLO and 
so forth). The customary fixed-order component [d(5‘aia2]fo consequently also the finite 

component) definitely contains azimuthal correlations with respect to although these are not 
logarithmically-enhanced in the small-gr region. 

To obtain NLL-I-NLO accuracy we have to include the functions and ^( 2 )(Afi,Af 2 ) 
generalized form factor of Eq. (0), the function 7 ^d)(^i.^ 2 ) in the hard/collinear factor 

finite component of Eq. (ITOll up to 0{as). To reach NNLL-I-NNLO 
accuracy we need to include also the functions g(3)(Afi,V2)^ finite component up 

to C>(q;|)@. This matching procedure between resummed and finite contributions guarantees to 
achieve uniform theoretical accuracy over the entire range of transverse momenta. In particular, 
we remark that the inclusion of resummed component at the NNLL-I-NNLO level 

is essential to achieve NNLO accuracy in the small-g^ region (considering a generic upper limit 
value gumax, the integral over the range 0 < < gTmax of the qt distribution at the NNLL-I-NNLO 

level includes the complete perturbative terms up to NNLO). An analogous remark applies to the 
inclusion of NLL-I-NLO level. 


We have so far illustrated the resummation formalism for the most general sixfold differential 


partonic cross section 


da 


aia 2 — 


(and for the corresponding hadronic cross section in Eq. (j3|)). 


Starting from 


da 


aia.2^^3/4 


and performing integrations over some kinematical variables, we can 


obtain resummed results for more inclusive g^-dependent distributions. For instance, integrating 
over the lepton kinematical variables fl, we obtain the qt cross section da/{d'^ qt dM'^ dy) at fixed 
invariant mass and rapidity of the lepton pair. The corresponding resummed component of the 
partonic cross section, as obtained from Eq. dSD, is 


^ aia2—>-/3/4 

qt dM'^ dy 


{Qt, M, y, s; as, t Wa,a2Mb2^v{qT, M, y, s; as, p^), 


bi,b2=qf,qfi 


( 11 ) 

where (T®q^,_^; 3 ^(M^) is the Born level (EW) total cross section for the partonic subprocess g/g/' 

V /3Z4. By performing an additional integration over the rapidity y of the vector boson (lepton 
pair), we obtain daj[d^q^dM'^) and the corresponding resummed component of the partonic cross 
section simply involves the integration over y of the resummed factor W{q^, M, y, s) in Eqs. ([5]) and 
dm) (or, equivalently, the factor W{h, M,y, s) in Eq. (|6])). After integration over y, the ensuing 


^This classification of the resummed-|-matched expansion exactly coincides with that of Refs. [26l|40]. We simply 
note that we are using labels that differ from those used in Refs. [26l|40]. The various terms of the expansion are 
denoted here (analogously to Ref. [43]) with the labels NLL-I-NLO and NNLL-I-NNLO, whereas they were denoted 
in Refs. [5S1|30] with the corresponding labels NLL-I-LO and NNLL-I-NLO. The fixed-order labels NLO and NNLO 
used here directly refer to the perturbative accuracy in the small-gr region (which corresponds to the perturbative 
accuracy of the total cross section), whereas the labels LO and NLO used in Refs. [26l [40] were directly referring 
to the perturbative accuracy in the large-gr region. 
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resummed factor depends on M and s, and it can be conveniently expressed in exponentiated 
form [26] by considering ‘single’ N Mellin moments with respect to the variable 2 ; = M'^/s at 
hxed M. The resummed expression for these ‘single’ N moments is exactly obtained by simply 
setting Ni = N 2 = N in Eqs. ©-(El). Our resummed calculation of da/{d? qt dM"^) was discussed 
in Ref. [39lll0|, and it is implemented in the numerical code DYqT. In Refs. [391110] we presented 
detailed quantitative results for vector boson production at Tevatron energies. Results from DYqT 
at LHC energies are presented in the following Sect. 13.11 

Within our formalism the resummation of the large terms ln(M^/g|,) at small values of qt 
is achieved by hrst performing the Fourier transformation of the qt cross section (or, more pre¬ 
cisely, of its singular behaviour in the small-g^r region) from space to b space (incidentally, the 
renormalization scale fin and the others auxiliary scales Q and fip are kept hxed and, especially, 
independent of qt in the integration over qt of the Fourier transformation). In b space, the large 
logarithmic variable (whose dependence has to be resummed) is L, at large values of b. Note that 
in the context of the resummation approach, the parameter as{fi‘n)L is formally considered to 
be of order unity. Therefore, the ratio of two successive terms in the expansion ([H]) is formally 
of O^asifiji)) (with no L enhancement). In this respect the resummed logarithmic expansion in 
Fq. (|H]) is as systematic as any customary hxed-order expansion in powers of asifiji)- Analo¬ 
gously to any perturbative expansions, the perturbative terms M/fin, M/Q) 

in Fq. ([8]) have an explicit logarithmic dependence on ln{M/fin) or ln(M/Q) (see, e.g., Fqs. (22) 
and (23) in Ref. [26]). Therefore, to avoid additional large logarithmic enhancements that would 
spoil the formal behaviour of the expansion in Fq. ([S]), the renormalization scale fin has to be set 
at a value of the order of M ~ Q. A completely analogous reasoning applies to the fip depen¬ 
dence of (M; asifiji), M/fin, M/fip, M/Q) in the expansion of Eq. ([9]) and, therefore, we 

should set Pi? ~ M. In other words, once the enhanced perturbative dependence on (i.e., 

on the two different scales M and 1/6) is explicitly resummed (albeit at a dehnite logarithmic 
accuracy), we are effectively dealing with a single-scale observable at the hard scale M and we can 
set pij ~ Pi? ~ M in both the resummed and hnite components of the qt cross section in Eq. (|1]). 

We remark that setting fi = 0{M) (here fi generically denotes the auxiliary scales fiji,fJ^F,Q) 
does not mean that the qt cross section is physically controlled by parton radiation with intensity 
that is proportional to The resummed form factor exp{^^^i’^^^} in Eq. © (and the 

ensuing logarithmic expansion in Eq. ©) is produced by multiparton radiation with intensity 
that is proportional to as{k‘^) and is a dynamical scale that varies in the range > k'^ > \jb^ 
(see, for instance, Eq. (19) in Ref. [2S]), where 1/6^ can be physically identihed with at small 
values of qn- Setting p ~ M in Eqs. © and © corresponds, roughly speaking, to consider the 
scale range fi^ > k‘^ > l/b‘^ (it does not correspond to set k'^ ^ f^ ^ M^). 

We recall [26] a feature of our resummation formalism. The small-gT singular contributions 
that are resummed in Eqs. © (or Eq. fllip l are controlled by the large logarithmic parameter 
ln(M^/g|,), which corresponds to L = ln{Q‘^b^/b^) (with Q ~ M) in b space at 6 —)■ cxo. In our 
resummation formula ©, we actually use the logarithmic parameter L = ln((5^6^/6g -|- 1) |41j . 
The motivations to use the logarithmic parameter L are detailed in Ref. [26] (see, in particular, 
the Appendix B and the comments that accompany Eqs. (16)-(18) and Eqs. (74)-(75) in Ref. |26j). 
and here we simply limit ourselves to recalling some aspects. In the relevant resummation region 
bQ 3> 1, we have L = L + 0{l/{Q‘^b‘^)) and, therefore, L and L are fully equivalent to arbitrary 
logarithmic accuracy (in other words, the replacement L L simply modihes the partition of 




small-gT non-singular contributions between the two components in the right-hand side of Eq. (P|)). 
However, L and L have a very different behaviour as 6 —)■ 0 (and, thus, they differently affect the 
qt cross section in the large-g^ regioiJ§). When 6Q <C 1, we have L S> 1 and, therefore, the re¬ 
placement L —)■ L in Eq. ([7]) would produce the resummation of large and unjustihed perturbative 
contributions in the large-g-r region (strictly speaking, the replacement L —)■ L leads to a qt cross 
section that is even not integrable over qt when qt —)■ oo: see, in particular, Eqs. (131) and (132) 
of the arXiv version of Ref. |26] and related accompanying comments). In contrast, when 6Q -C 1 
we have L —?• 0 and ^(^ 1 X 2 ) q Therefore, the use of L reduces the impact of unjustified large 
contributions that can be introduced in the small-6 region through the resummation procedure. 
Moreover, the behaviour of the form factor at 6 = 0 is related to the integral over qt 

of the g^-dependent cross section and, since we have exp{^(^i’'^ 2 )| _ resummation 

formalism fulhls a perturbative unitarity constraint |26]: after inclusion of the hnite component 
as in Eq. m. the integration over qx of our resummed qt cross sections recovers the hxed-order 
predictions for the total cross sections. Specihcally, the integral over qx of da/ {d? qt dM‘^dy) and 
da/{d?Q t dM"^) at the NNLL-I-NNLO (NLL-I-NLO) accuracy completely and exactly (i.e., with 
no additional higher-order contributions) agrees with the rapidity distribution da/{dM‘^dy) and 
the total cross section da/dM'^ at NNLO (NLO) accuracy, respectively. In summary, the expres¬ 
sions ([7]) and dH]) in terms of the logarithmic parameter L correctly resum the large parametric 
dependence on ln(6Q) at large values of bQ and they introduce parametrically-small perturbative 
contributions at intermediate or small values of bQ (the coefficients of the perturbative corrections 
are proportional to powers of L with L (9(1) if bQ ~ (9(1) or L 1 if 6Q 1). After having 
combined the resummed calculation at N^LL accuracy with the complete N^LO calculation, as 
in Eqs. (jl|) and flTOl) . these parametrically-small corrections produce residual terms that start to 
contribute at the level. Therefore, the use of L has the purpose of reducing the impact of 

unjustihed and large higher-order (i.e., beyond the N^LO level) contributions that can be possibly 
introduced at intermediate and large values of qt through the resummation of the logarithmic 
perturbative behaviour at small values of qt- In particular, no residual higher-order contributions 
are introduced in the case of the total (integrated over qx) cross section (which is the most basic 
quantity that is not affected by logarithmically-enhanced perturbative corrections). 

We add some relevant comments about the dependence of the resummed cross section on the 
kinematical variables that specify the angular distribution of the leptons with respect to the 
vector boson. By direct inspection of Eqs. ([S]) and fllll) we see that they involve exactly the same 
resummation factor W. The only difference between the right-hand side of these equations arises 
form the Born level factors dd^^'* /dO, and which are related as follows through the integration 

over ri: 

dd^^^ 

M\ il) , ( 12 ) 

with the normalization condition 

j dfl Fq^q^,^i^i^{qx/M; f2) = 1 . (13) 

Although both factors depend on EW parameters (EW couplings, mass and width of the vector 
boson), they have a different dependence on the relevant kinematical variables. The vector boson 

§The contribution of the integral in Eq. ([6]) from the integration region where h<0{l/M) ^ 0{1/Q) always 
gives (provided yV(b,M) is integrable over such region) a non-singular contribution to the qx cross section in the 
small-gr region. 
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distribution {(PqxdM"^ dy) (and, analogously, /{d?qxdM"^)) involves the Born level 

total cross section which depends on M^, whereas the less inclusive leptonic distribution 

dy involves the Born level differential cross section da^^^dfl that addition¬ 
ally depends on f 2 and also on the transverse momentum qx of the lepton pair (see the function 
F in the right-hand side of Eq. (IT^ ). To our knowledge the qx dependence of da^^l/dCt has not 
received much attention in the previous literature on transverse-momentum resummation and, 
therefore, we discuss this issue with some details in Appendix Physically, this qx dependence 
is a necessary consequence of transverse-momentum conservation and it arises as a gr-recoil effect 
in transverse-momentum resummation. At the LO in perturbation theory the lepton angular dis¬ 
tribution is determined by the Born level production and decay process of the vector boson, which 
carries a vanishing transverse momentum. Through the resummation procedure at fixed lepton 
momenta, higher-order contributions due to soft and collinear multiparton radiation dynamically 
produce a finite value of the transverse momentum qx of the lepton pair, and this finite value of 
qx has to be distributed between the two lepton momenta by affecting the lepton angular distri¬ 
bution. This gr-recoil effect on the Born level angular distribution is a non-singular contribution 
to the qx cross section at small values of qx and, therefore, it is not directly and unambiguously 
computable through transverse-momentum resummation. In other words, the Born level function 
F in Eq. (IT^ has the form 

Fqfq^i^i^uiqx/M] M^, f2) = M^, fl) + 0{qx/M) , (14) 

where E(0; M^, il) is uniquely determined, whereas the small-gr corrections of 0{qx/M) has to be 
properly specified. In any physical computations of lepton observables (i.e., in any computations 
that avoid possible unphysical behaviour due to violation of momentum conservation for the 
decay process q = p^{h) + through transverse-momentum resummation, a consistent qx- 

recoil prescription has to be actually (either explicitly or implicitly) implementecfl. Note that, 
after having combined the resummed and finite components as in Eqs. (jl]) and (ITOD . the 0{qx/M) 
recoil effects lead to contributions that start at 0(a|) (i.e., N^LO) in the case of our resummed 
calculation at NNLL-I-NNLO accuracy (correspondingly, these contributions start at 0(a|) in the 
case of NLL-I-NLO accuracy). Obviously there are infinite ways of implementing the ^-recoil 
effect, and in Appendix we explicitly describe a very general and consistent procedurqllj. Note 
that the gr-recoil effect completely cancels after integration over the leptonic variables (see 
Eq. dH). 

Our resummed calculation of the sixfold differential distribution in Eq. ([2]) is implemented in 
the numerical code DYRes, which allows the user to apply arbitrary kinematical cuts on the mo¬ 
menta of the final-state leptons and to compute the corresponding relevant distributions in form of 
bin histograms. We add some comments on the numerical implementation of our calculation. In 
Eqs. 0-0 we have illustrated the structure of the resummed component in the double (A^i, N 2 ) 

^The dynamical treatment of the qx recoil is embedded in the formulation of transverse-momentum (kx) fac¬ 
torization [ 8 OI [81] of hard-scattering processes at high energy (at small x). 

II The QT-recoil issue is not a specific issue of leptonic decay in vector boson production. The issue is completely 
general (see Appendix [^ , and it arises in any qx resummed calculation for the production of a set of particles 
with measured momenta at fixed total transverse momentum qx (e.g., diphoton, diboson, or heavy-quark pair 
production). A noticeable exception (as discussed in Appendix lA)) is the production of a SM Higgs boson and its 
subsequent decay. In this case, due to the spin-0 nature of the Higgs boson, the qx dependence of the correspond¬ 
ing Born level function F{qr/M; M'^,n) can be entirely determined by kinematics [13], without the necessity of 
specifying gT’-recoil effects of dynamical origin. 
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Mellin space. Through the inverse Mellin transformation, this structure can equivalently be ex¬ 
pressed in terms of convolutions with respect to longitudinal momentum fractions Xi and X 2 (see 
Eq. ([3])). In the DYRes code, the Mellin inversion is carried out numerically. The results for the 
NNLO term 'Hy'’ in Eq. (|9]) are presented in Ref. [71] in analytic form directly in {xi,X 2 ) space. 
These results have to be transformed in Mellin space. Then, the Mellin inversion requires the nu¬ 
merical evaluation of some basic iV-moment functions that appear in the expression of 
this evaluation has to be performed for complex values of iV, and we use the numerical results of 
Ref. [82|. This implementation of the resummed component is completely analogous to that of the 
DYqT code |39l HQ] and of other previous computations Nonetheless, the efficient generation 
of ‘vector boson events’ according to the multidifferential distribution of Eq. (jl]) and the inclusion 
of the leptonic decay are technically non trivial, and this requires substantial improvements in the 
computational speed of the numerical code that evaluates the resummed component of the cross 
section. The hxed-order (NLO and NNLO) cross section in Eq. flTOD and then the hnite component 
of the cross section in Eq. (|3]) are evaluated through an appropriate modihcation of the DYNNLO 
code [5]: DYNNLO is particularly suitable to this purpose, since it is based on the qt subtraction 
formalism [23], which uses the transverse-momentum resummation formalism to construct the 
subtraction counterterms. 


Using the resummation expansion parameter L in Eq. (|7]) and the matching procedure (which 
implements the perturbative unitarity constraint on the total cross section) with the complete 
hxed-order calculation, our resummation formalism [26] formally achieves a uniform theoretical 
accuracy in the region of small and intermediate values of qt, and it avoids the introduction of 
large unjustihed higher-order contributions in the large-gr region. In the large-gr region, the 
results of the resummed calculation are consistent with the customary hxed-order results and, 
typically [26l HQ], show larger theoretical uncertainties (e.g., larger dependence with respect to 
auxiliary-scale variations) with respect to the corresponding hxed-order results. This feature is not 
unexpected, since the theoretical knowledge (and the ensuing resummation) of large logarithmic 
contributions at small qt cannot improve the theoretical predictions at large values of qx- In the 
large-gy region, where the resummed calculation shows ‘unjustihed’ large uncertainties and ensuing 
loss of predictivity with respect to the hxed-order calculation, the reliability of the resummed 
calculation is superseded by that of the hxed-order calculation. In this large-gr region, we can 
simply use the theoretical results of the hxed-order calculation. In the computation of quantities 
that directly and explicitly depend on qt (e.g., the transverse-momentum spectrum of the vector 
boson), it is relatively straightforward to identify and select ‘a posteriori’ the large-g^ region 
where the resummed calculation is superseded by the hxed-order calculation. In the present work, 
however, we are also interested in studying kinematical distributions of the vector boson decay 
products: our goal is thus to generate the full kinematics of the vector boson and its (leptonic) 
decay, to apply the required acceptance cuts, and to compute the relevant distributions of the 
lepton kinematical variables. In this framework, the actual results can become sensitive to the 
large-gr region in which the resummed calculation cannot improve the accuracy of the hxed-order 
calculation. To reduce this sensitivity, in the DYRes implementation of the resummed calculation 
we thus introduce a smooth switching procedure at large value of qt by replacing the resummed 
cross section in Eq. (|H]) as follows: 


da, 


aia.2~^^3U 


-)■ w(gr) 


dd 


(res.) 

aia2—>-/3/4 


+ 


dd 


(fin.) 

aia2^l3l4 


j -L (1 - w{qT)) dd, 


^aia2—>-/3U 


J f.o. 


( 15 ) 
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where the function w^qx) is defined as 

( 16 ) 


( 17 ) 

We have quantitatively checked that the value of the parameter q^' can be ‘suitably’ chosen in 
the large-gT region, and that both parameters q^' and Aqx can be consistently chosen so as not to 
spoil our unitarity constraint (in Sect. 13.II we show that the integral over qx of our NLL+NLO and 
NNLL+NNLO resummed results still reproduces well the NLO and NNLO total cross sections). 
We note that we do not introduce any switching procedure in the DYqT calculation (though, its 
introduction is feasible) since, as previously mentioned, the identihcation of the large-g^ region is 
straightforward in the computation of da/{d?qx dM"^). 

We recall [26] that the resummed form factor exp{^(a 5 ,L)} of Eq. (I7|) is singular at very 
large values of b. The singularity occurs in the region b> 1/Aqcd, where Aqcd is the momentum 
scale of the Landau pole of the perturbative running coupling ^^(/i^). This singularity is the 
‘perturbative’ signal of the onset of non-perturbative (NP) phenomena at very large values of b 
(which practically affect the region of very small transverse momenta). In this region NP effects 
cannot any longer be regarded as small quantitative corrections and they have to be taken into 
account in QCD calculations. A simple and customary procedure to include NP effects is as 
follows. The singular behaviour of the perturbative form factor exp{^(Q; 5 , L)} is removed by 
using a regularization procedur^ and the resummed expression in Eq. ([7|) is then multiplied by a 
NP form factor and it is inserted as integrand of the b space integral in Eq. (Ej). The regularization 
procedure that was used in the DYqT calculation |10| is the ‘minimal prescription’ of Ref. [SH l6T] . 
which basically amounts to avoid the singularity of exp{^(Q; 5 , L)} by deforming the integration 
contour of Eq. ([6|) in the complex b plane. In the DYRes calculation of the present work, we use a 
different regularization procedure by freezing the b dependence of exp{^(a 5 ', L)} before reaching 
its singular point. The freezing procedure follows the ‘6* prescription’ of Refs. [221 123 ] and it is 
obtained by performing the replacement 

b^^bl = b^ blJib^ + blJ (18) 

in the b dependence of Q{as,L). The value of the parameter 6iim has to be large (&iimAf ~ 
b\imQ 3> 1) but smaller than the value of b at which the singularity of exp{^(Q( 5 , L)} takes place 
(note that the replacement in Eq. ([T8|) has a negligible effect at small and intermediate values 
of b since bl = 6^(1 + 0{b‘^Q'^/b^^^Q'^)) ~ 6^ if bQ < 1). The use of the 6* freezing procedure 
improves the (numerical) performances of the DYRes code. Additional comments on NP effects 
are presented in Sect. 13.11 

**We recall that the resummed form factor exp{C/(as, Y)} produces a strong suppression (Q{as,L) oc —asLF') 
in the large-6 region where asLd' ^ (D(l). Therefore, the choice of different regularization procedures mildly affects 
[inillllllS] the results since its effects are relevant only in the region 6 ^ 0(1/ Aqcd) where the 6 integral is strongly 
damped by the form factor. 


w{qx) = 

and the function f{qx) is chosen as 


qr < q^- 


fiqx) qr > qx 


fiqx) = exp - 


(g^- - QtY 

2{Aqxy 
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3 Numerical results at the LHC 


In this Section we consider the processes pp —)■ Zj^* —)■ and pp —)■ —)■ Ivi at LHC 

energies. We present our resummed results at NNLL+NNLO and NLL+NLO accuracy, and we 
compare them with some of the available LHC data. We compute the hadronic cross sections at 
NNLL+NNLO (NLL+NLO) accuracy by using the NNPDF3.0 NNLO (NLO) [HS] parton densities 
functions (PDFs), with as{m\) = 0.118 and with as{p%) evaluated at 3-loop (2-loop) order. As 
in customary hxed-order calculations at high invariant mass (M = 0{mz)), we consider Nf = 5 
flavours of light quarks and we treat them in the massless approximation. 

As for the EW couplings, we use the so called scheme, where the input parameters are 
Gf, rnz, mw- In particular, we use the PDG 2014 [HS] values Gp = 1.1663787 x 10“^ GeV“^, 
mz = 91.1876 GeV, Tz = 2.4952 GeV, mw = 80.385 GeV and Tw = 2.085 GeV and in the 
case of production, we use the (unitarity constrained) GKM matrix elements Vud = 0.97427, 
Vus = 0.22536, 14b = 0.00355, 14d = 0.22522, 14* = 0.97343, 14b = 0.0414. Our calculation 
implements the leptonic decays Z/y* —)• and W Ivi (we include the effects of the 
interference and of the hnite width of the W and Z bosons) with the corresponding spin correlations 
and the full dependence on the kinematical variables of hnal state leptons. This allows us to 
take into account the typical kinematical cuts on hnal state leptons that are considered in the 
experimental analysis. As discussed in Sect. [2], the resummed calculation at hxed lepton momenta 
requires a gy-recoil procedure. We implement a procedure that is described in Appendix and 
that is practically equivalent to compute the Born level distribution /dfi of Eqs. ([5]) and flT^ 
in the Gollins-Soper rest frame [TD] (this is exactly the same procedure as used in other resummed 
calculations [551 [561 EH [52]). As explained in Sect. [2l the DYRes resummed calculation uses a 
smooth switching procedure (see Eq. flT^ l in the large-gr region. In our numerical implementation 
the parameters in Eq. f[T7|) are chosen to be Agx = M/(2\/2) and = 3M/4. 

3.1 Inclusive results at fixed qp 

We start the presentation of our results by discussing some general features of the gp spectrum that 
can be addressed at the inclusive level, i.e. after integration over the lepton angular variables and 
over the rapidity y of the lepton pair. Unless otherwise specihed, the numerical results of this Sub¬ 
section are obtained by using the code DYqT [321 SO]- The code DYqT is publicly available and it can 
be downloaded from the URL address http://pcteserver.nii.infn.it/~ferrera/dyqt.html . 

We hrst consider the dependence on the auxiliary scales gip, and Q. These scales have to 
be set at values of the order of the invariant mass M of the produced system, with no dehnite 
preference for specihc values. Then, scale variations around the chosen central value can be used to 
estimate the size of yet uncalculated higher-order terms and the ensuing perturbative uncertainties. 
In the NNLL+NNLO studies of Refs. [26102103] on Higgs boson production, and in our previous 
work on vector boson production HI the central reference values of the scales were chosen as 
Pf = Pr = 2Q = mpI where is the mass of the produced boson (the Higgs boson mass in the 
case of Higgs boson production, and the vector boson mass in the case of vector boson production). 
In the case of Higgs boson production, this choice gives consistent NLL+NLO and NNLL+NNLO 
results with a reduced scale dependence at NNLL+NNLO level and, in particular, with a nice 
overlap of the NLL+NLO and NNLL+NNLO uncertainty bands (see, e.g.. Fig. 2 of Ref. [12])- In 
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RATIO da/dq^ (pb/GeV) 


the case of vector boson production, our previous studies were focused on Tevatron energies, and 
a similar pattern was observed [iQ]. When moving to vector boson production at LHC energies, 
we notice that the factorization-scale dependence exhibits a (slightly) different behaviour. 




qT (GeV) qT (GeV) 

(a) (b) 

Figure 1: The qt spectrum of Z bosons at the LHC (y/s = 8 TeV). The bands of the NLL+NLO 
(red dashed) and NNLL+NNLO (blue solid) results are obtained by performing fip variations 
around the central values fip = (l^ft panel) and (right panel). The lower panel 

presents the ratio of the NLL+NLO and NNLL+NNLO results with respect to the NNLL+NNLO 
result at the central value of pp. The result of the convolution of the NLL+NLO partonic cross 
section with NNLO PDFs at the central value of pp is also reported (black dotted) in the lower 
panel. 

In Fig. [T] we present results for the qp spectrum of on-shell Z bosons produced at the LHC 
[y/s = 8 TeV) and the corresponding dependence on factorization-scale variations. In the left 
panel the central scale is pp = pp = 2Q = mz, while in the right panel the central scale is 
Pf = Pr = Q = "nizl^. In both panels, we present the NLL-I-NLO and NNLL-I-NNLO results 
at the central scale and corresponding bands that are obtained by varying (up and down) the 
factorization scale by a factor of 2 around its corresponding central value. The lower panel in 
Figs. [11(a) and [11(b) presents the ratio of the various results with respect to the NNLL-I-NNLO 
result at the corresponding central scale. \i pp = pp = 2Q = mz is the central scale choice (left 
panel), we see that the factorization-scale bands at NLL-I-NLO and NNLL-I-NNLO accuracy never 
overlap, except for the tiny region around qp ^ 7 GeV where they cross each other. The lack of 
overlap is particularly evident in the peak region, where the bulk of the events is placed: here the 
central NLL-I-NLO and NNLL-I-NNLO results differ by about 30%. We also notice that throughout 
the region of small and intermediate values of qp {qp ^ 30 GeV) the size of the NLL-I-NLO band 
is rather small and it is always (with the exception of a small region around qp ^ 8 GeV) smaller 
than the size of the NNLL-I-NNLO band, and this suggests that an accidental cancellation of 
the Pp dependence may occur at the NLL-I-NLO level with this choice of central scale. In the 
right panel we observe a /Xir-dependence behaviour that is qualitatively similar but quantitatively 
different from that in the left panel, li pp = pp = Q = mzl 2 is the central scale choice (right 
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panel), the NLL+NLO and NNLL+NNLO bands are closer and they overlap at small transverse 
momenta. The overlap occnrs in a limited region of qt that, nevertheless, inclndes the peak region. 
The shape of the spectra appears closer when going from NLL+NLO to NNLL+NNLO accnracy, 
and the NLL+NLO band is wider than the NNLL+NNLO one in the small and intermediate 
region of qt- Note that the central valnes of fiR are fiR = irtz and ^r = irtzl^ in the left and right 
plot, respectively, bnt we have checked that this difference has little effect: the observed different 
behavionr is mainly dne to the different central valne of In snmmary, the /ij? dependence 
observed in the left panel of Fig. [1] snggests that the corresponding scale variation bands (and 
especially the NLL+NLO band) are likely to nnderestimate the pertnrbative nncertainties of the 
calcnlation. Based on these observations, in the rest of the paper, we will adopt = fJ^R = Q = 
mzl‘2^ as reference valnes of the central scales. 

The differences between the NLL+NLO and NNLL+NNLO resnlts have a twofold origin. Part 
of the differences is dne to the next-order radiative corrections in the partonic cross sections, and 
the remaining part is dne to the increased order of the PDFs. To qnantify the impact of these 
two different contribntions, we have considered the result that is obtained by convoluting the 
NLL+NLO partonic cross sections with the NNLO PDFs. This result, at the central scales that 
are considered in Figs. [I](a) and [11(b), is reported (see the black dotted line) in the corresponding 
lower panel, and we can see that it is quite close to the NLL+NLO result with NLO PDFs. In 
other words, a large part of the quantitative differences between the NLL+NLO and NNLL+NNLO 
results is due to the corresponding differences at the level of the partonic cross sections. 




qT (GeV) qT (GeV) 

(a) (b) 

Figure 2: NLL+NLO (red dashed) and NNLL+NNLO (blue solid) results for the qt speetrum of 
Z bosons at the LHC with energies s/s =8 TeV (left panel) and ^/s =14 TeV (right panel). The 
NLL+NLO and NNLL+NNLO results with the central value fip = P>r = Q = ''tiz/2 of the scales 
are enclosed by corresponding bands. The bands are obtained by performing {fip, Q} variations 
(as described in the text) around the central value razf^. The lower panel presents the ratio of the 
scale-dependent NLL+NLO and NNLL+NNLO results with respect to the NNLL+NNLO result at 
the central value of the scales. 
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The NLL+NLO and NNLL+NNLO results for the qx spectrum of on-shell Z boson produced at 
the LHC with different collision energies are presented in Fig. |2j We consider two centre-of-mass 
energies: y/s =8 TeV (Fig. |2]left) and \fs =14 TeV (Fig. [2] right). At each logarithmic accuracy 
we present the result at the central value fix = = Q = 'mzl‘2. of the scales and a corresponding 

band. The bands provide an estimate of the perturbative uncertainties of the calculations due 
to missing higher-order contributions. The bands are obtained through independent variations of 
Pi?, pij and Q by following the procedure of Ref. [iQ]: we independently vary pi?, pij and Q in 
the range mz!^ < {pi?,pii, Q} < mz with the constraints 0.5 < < 2 and 0.5 < Q/iir < 2. 

We remind the reader that the constraint on [ixlii-R is introduced to avoid large logarithmic 
contributions (ln(pi?/pij) terms from the evolution of the parton densities) in the perturbative 
expansion of the hard/collinear factor "Hy of Eq. (JTI). Analogously, the constraint on Q/pi^ avoids 
large logarithmic terms (ln(Q/pii)) in the resummed expansion of the form factor exp{^} of 
Eq. d?]). The lower panels in Fig. [2] present the ratio of the scale-dependent NLL-I-NLO and 
NNLL-I-NNLO results with respect to the NNLL-I-NNLO result at the central value Pi? = p_R = 
Q = mzl‘2^ of the scales. 

The region of small and intermediate values of qt is shown in the main panels of Fig.O At hxed 
centre-of-mass energy the NNLL-I-NNLO qx spectrum is harder than the spectrum at NLL-I-NLO 
accuracy. At fixed value of qx the cross section sizeably increases by increasing the centre-of-mass 
energy from 8 TeV to 14 TeV. The shape of the NNLL-I-NNLO qx spectrum is slightly harder at the 
higher energy. The NLL-I-NLO scale-variation band is wider than the NNLL-I-NNLO band. The 
NLL-I-NLO and NNLL-I-NNLO bands overlap at small transverse momenta and remain very close 
by increasing qx (the differences with respect to the plot on the right-hand side of Fig. [T]are due 
to the additional dependence on p/j and Q). The NNLL-I-NNLO (NLL-I-NLO) scale dependence 
is about ±10% (±20%) at the peak, it decreases to about ±2% (±7%) at qx — 10 GeV and 
increases to about ±6% (±10%) at qx ~ 25 GeV. Since the NNLL±NNLO and NLL±NLO bands 
do not exactly touch each other in the region where gr ~ 8 GeV, one may argue that the ‘true’ 
perturbative uncertainty of the NNLL±NNLO result in this region is slightly larger than the size of 
the NNLL±NNLO scale dependence band (for instance, one may use [39] the difference between 
the NNLL±NNLO central scale result and the upper line of the NLL±NLO band in Fig. [2] to 
estimate the uncertainty of the NNLL±NNLO result). 

The inset plots show the cross section in the large-g^ region. The resummation results ob¬ 
tained with DYqT and reported in the inset plots are presented for completeness and mainly for 
illustrative purposes. At large values of qx {qx ^'^z) the resummed result looses predictivity, and 
its perturbative uncertainty becomes large. In this region of transverse momenta we see that the 
uncertainty band increases in going from the NLL±NLO to the NNLL±NNLO level. However, 
as already mentioned in Sect. [21 at high qx the resummation cannot improve the predictivity of 
hxed-order calculations and the DYqT result in Fig. [2] cannot be regarded as reference theoretical 
result. The resummed result has to be replaced by the standard hxed-order prediction. The 
NNLO (NLO) result (which is not shown in Fig. [2]) lies inside the NNLL±NNLO (NLL±NLO) 
band and the former has a smaller scale dependence than the latter. We also note that, at high 
qx, the preferred reference central scales hr and pi? of the hxed-order prediction should be of the 
order of (rather than of the order of mz)- 

We also recall that, increasing qx throughout the high-g^ region, hxed-order QGD calculations 
are ahected by additional and potentially-large logarithmic terms. These are collinear (fragmen¬ 
tation) contributions [87l[88], which become more relevant by increasing the ratio qx/M, and soft 
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(threshold) contributions [5911^ . which become more relevant by increasing the ratio qr/y/s (or 
Qt/'/S). 

We have so far discussed only uncertainties from missing higher-order contributions. Before 
moving to consider the case in which cuts on the hnal-state leptons are applied, we briefly discuss 
two additional sources of QCD uncertainties on the qt spectrum: the uncertainty from PDFs and 
that from NP effects. We consider these effects in turn. 

Modern sets of PDFs include an estimate of the errors (mainly experimental errors) in their 
determination from global data hts, and this estimate can then be used to compute the ensuing 
PDF uncertainty on the QCD calculation of hadron collider observables. In Fig. [3] we consider Z 
boson production at NNLL-I-NNLO accuracy. In Fig. [3] (a) we report the NNLL-I-NNLO results of 
Fig. [5] (a) (v^ =14 TeV) and the effect of the PDF uncertainty at 68% CL on the NNLL-I-NNLO 
calculation at the central scale value fip = fJ-R = Q = tuz/^. In Fig. |3](b) the scale-dependence and 
PDF-uncertainty bands are normalized to the central NNLL-I-NNLO prediction, and we present 
results at both energies y/s =8 TeV (lower panel) and yfs =14 TeV (upper panel). We see that 
the PDF uncertainty is smaller than the scale uncertainty. Moreover, the PDF uncertainty is 
approximately independent on the transverse momentum, and it has a value of about ±3% at 
both energies y/s = 8 TeV and y/s = 14 TeV. 

NP effects are known to increasingly affect the transverse-momentum spectrum as qp decreases 
towards g-r —^ 0. A detailed study of these effects is beyond the scope of the present work. We 
limit ourselves to roughly estimate the possible impact of such effects, and we use a very simple 
model in which the perturbative form factor exp{^(a 5 ',L) in Eq. ([7]) is multiplied by a NP form 
factor SNp{b) = exp{—g]yfpb^}, which produces a Gaussian smearing of the qp distribution at 
small-gT values. We vary the value of the parameter g^p in a quite wide (‘conservative’) range, 
0 < gNP < 1-2 GeV^, and in Fig. [3](a) (black band) we show the ensuing quantitative effects 
on the qp spectrum. In Fig. [3](b) the NP effects are normalized with respect to the perturbative 
NNLL-I-NNLO result at central value of the scales. 

Gomparing the lower panels of Fig. [2] with Fig. [3], we can hrst make an overall qualitative 
comment. Perturbative corrections make the qp spectrum harder in going from NLL-I-NLO to 
NNLL-I-NNLO accuracy, and this occurs at both small and intermediate values of qp. NP effects 
increase the hardness of the qp spectrum at small values of qp and they are negligible at interme¬ 
diate values of qp. Therefore, we note a non trivial interplay of perturbative and NP effects. In 
particular, at small values of qp higher-order perturbative contributions can be mimicked by NP 
effects. 

At the quantitative level, in Fig. [3] we see that the NNLL-I-NNLO result supplemented with 
NP effects is very close to the perturbative result except in the very low qp region {qp^3 GeV), 
i.e. below the peak of the qp distribution. In the region 3 GeV gr ~ 10 GeV, the size of the NP 
band is similar to that of the PDF uncertainty band. At larger values of qp, the NP effects vanish 
(the size of the NP band is smaller than about 2% starting from qp ~ 15 GeV). 

We note that our simple model treats the regularization of the perturbative form factor 
(through the ‘minimal prescription’, see Sect. [2]) and the NP form factor in an uncorrelated 
way, and this produces a conservative estimate of NP uncertainties. In other words, the model 
underestimates the potential of the resummed calculation at very small values of qp. For instance. 


17 


the NP model can be improved by correlating the interplay between the perturbative form factor 
(and, e.g., its scale variation dependence) and the NP form factor (and the value of gNp), and 
further constraints on the NP model can be possibly obtained by inputs from comparisons with 
experimental data. 

In summary, from our brief discussion on the possible impact of NP effects for vector boson 
production at the LHC, we conclude that our conservative estimate leads to quantitative effects 
that are small and well within the scale variation dependence, still in the very low qt region. 
A quantitatively similar conclusion applies to the effect of PDF uncertainties. Based on these 
observations (and for practical purposes), in the presentation of our results of Sect. 13.21 we limit 
ourselves to considering only the perturbative calculation and the corresponding scale variation 
uncertainties. 



qi (GeV) 
(a) 



5 10 15 20 25 30 


qT (GeV) 

(b) 


Figure 3: (a) The qt spectrum at NNLL+NNLO accuracy for Z boson production at the LHC 
with y/s =14 TeV. Comparison of scale dependence (blue solid) and PDF (red crossed solid) 
uncertainties. The possible impact of NP effects is also shown (black crossed dashed), (b) The 
same results are normalized to the central NNLL+NNLO prediction at ^/s = 14 TeV (upper panel), 
and corresponding results are shown at y/s = 8 TeV (lower panel). 


We conclude this subsection by presenting a comparison between the DYqT results and the 
results of the ‘multidifferentiah program DYRes. When no cuts are applied on the final-state 
leptons, the qt spectrum of the on-shell vector boson obtained with DYRes has to be in agreement 
with the one obtained with the numerical program DYqT. We have numerically checked that this 
is indeed the case. For illustrative purposes, we show the results of a comparison in Fig. 01 Here 
we consider the qt spectrum for on-shell Z boson production at the LHC with ^/s = 7 TeV. The 
DYqT (solid line) and DYRes (histogram) results at central value of the scales are compared at 
both NLL-I-NLO (red) and NNLL-I-NNLO (blue) accuracy. At small and intermediate values of 
Qt (main plot in Fig. 0]), the DYqT and DYRes results agree (within the statistical uncertainties of 
the DYRes cod^) at both level of logarithmic accuracy. The quantitative degree of agreement is 

^^Here and in the following the errors reported in the tables and on the histograms refer to a numerical estimate 
of the accuracy of the Monte Carlo integration in the DYNNLD and DYRes codes. 
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more clearly visible in the lower panel, which presents the result of the calculation of the binned 
ratio between the DYRes and DYqT results at both NLL+NLO and NNLL+NNLO accuracy. The 
ratio is everywhere consistent with unity within the numerical accuracy of its computation (the 
numerical errors in the computation of the binned ratio are below about 1% at small values of 
qt, and they are still below about 2% in the region 30 GeV ^ gr ~ 50 GeV where the value of the 
cross section sizeably decreases). 



0 10 20 30 40 50 

Qt (GeV) 


Figure 4: The qt spectrum of on-shell Z bosons at the LHC (y/s = 7 TeV). Comparison 
of the DYRes (histograms) and DYqT (solid lines) resummed results at NLL-hNLO (red) and 
NNLL-hNNLO (blue) accuracy. The corresponding fixed-order calculations (dashed lines) at 0{as) 
(red) and C>(a|) (blue) are also shown. The ratio between the DYRes and DYqT results is shown in 
the lower panel. 

We recall (see Sect. E]) that, at the inclusive level, the DYqT and DYRes calculations involve 
differences in the numerical implementation and two additional differences related to the treat¬ 
ment of the very low qt and high-grp regions. At very low values of qt, the difference is due to 
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the regularization procedure of the perturbative form factor for very large values of the impact 
parameter b: the DYqT calculation uses the ‘minimal prescription’, while the DYRes calculation uses 
the 6* freezing procedure. In our actual calculation with DYRes the value of 6iim in Eq. fflSj) is set 
to 6iim = femax, wliere 6max E the maximum value of b that can be reached before encountering the 
singularity of the perturbative form factor (setting feiim = femax we do not introduce any additional 
regularization parameter, analogously to the case of the ‘minimal prescription’). The value of 6max 
depends on the renormalization and resummation scales fiR and Q and, in the case of Z and W 
production around the central value of the scales, the typical value is 6maxQ ~ 1-2 ■ We 

have checked that the ‘minimal prescription’ and the choice foum = &max give basically the same 
numerical results, also at very small values of qx {qx ~ 1 GeV). This numerical agreement is also 
visible (lower panel in Fig. 01) from the ratio between the DYRes and DYqT results at low values of 
qx- 


At large values of qx, the DYRes calculation implements the smooth switching procedure of 
Eqs. (IT^ - flTTI) . The large-gr region is shown in the inset plot of Fig. 01 and here the differences 
between the DYqT and DYRes calculations are due to the smooth switching procedure. In the 
high-gr region the DYRes result at NNLL+NNLO (NLL+NLO) accuracy basically agrees with 
the customary hxed-order result at 0(a‘^) (0{as))- The differences between the DYRes and DYqT 
results (consistently) decrease in going from NLL+NLO to NNLL+NNLO accuracy, and they are 
small at the NNLL+NNLO level. At both level of logarithmic accuracy, the DYRes and DYqT results 
agree within their corresponding scale variation uncertainties (which are not shown in the inset 
plot), and the DYRes result has a reduced scale dependence (it matches the scale dependence of 
the corresponding hxed-order result). The introduction of the smooth switching procedure in the 
DYRes calculation has practically a negligible quantitative effect on the unitarity constraint that 
is fulhlled by the DYqT calculation. In Tabled] we report the total cross sections for both Z and W 
production at ^/s = 7 TeV, and we compare the resummed DYRes results with the corresponding 
hxed-order results obtained with the DYNNLO code. We see that the NLL+NLO (NNLL+NNLO) 
total cross section agrees with the NLO (NNLO) result to better than 1% accuracy. 


Gross section [pb] 

NLO 

NLL+NLO 

NNLO 

NNLL+NNLO 

pp ^ Z ^ 

904.3 + 0.2 

904.6 ± 0.4 

949.1 + 0.7 

947.3 + 0.9 

pp —)■ 1+^^^ —)■ 

9819 ± 2 

9813 + 4 

10337 + 6 

10328 + 9 


Table 1: Total cross sections at the LHC (^/s = 7 TeV): fixed-order results and corresponding 
resummation results of the DYRes numerical program. 

In the main plot of Fig. 0] we also present a complementary information on the results of 
the hxed-order calculations (dashed lines) at 0{as) (red dashed) and C>(a|) (blue dashed). At 
intermediate values of qx the diherences between the resummed results at two subsequent orders 
are smaller than the diherences between the corresponding hxed-order results at two subsequent 
orders. The diherences between the resummed results and the corresponding hxed-order results 
sizeably increase by decreasing qx- At small values of qx, the result at 0{as) increases towards 
large positive values (they are outside the vertical size of the plot) and, in a hrst very small bin 
(not shown in the plot) around qx = 0, the 0{as) result would be very large and negative. The 
result at 0 {ag) has a very high unphysical peak (it is outside the vertical size of the plot) around 
qx 4 GeV, then it decreases towards very large negative values and, in a hrst very small bin 
(not shown in the plot) around qx = 0, the C>(a|) result would be very large and positive. 
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3.2 Vector boson production at the LHC 

In this Section we consider related) physical observables that depend on the individual lepton 
momenta and on the kinematics of the lepton pair. The dependence can be indirect, through the 
application of acceptance cuts, and direct, through the dehnition of the observable. Therefore, 
the resummed calculation presented in this Section are performed by using the numerical program 
DYRes. 

We start our presentation by considering the measurements of the qt spectrum of dilepton pairs 
at the LHC with ^/s = 7 TeV, as reported by the CMS [91] and ATLAS [92] Collaborations with 
an integrated luminosity of 36 pb“^ and 4.7 fb“^, respectively. The cuts that dehne the fiducial 
region in which the measurements are performed (our corresponding resummed calculation of the 
Z/y* spectrum is carried out in the same region) are as follows. In the case of the CMS analysis 
the invariant mass mu of the lepton pair is required to be in the range 60 GeV < mu < 120 GeV, 
and the leptons must be in the central rapidity region, with pseudorapidity \rf‘\ < 2.1, and they 
have a transverse momentum > 20 GeV. In the case of the ATLAS analysis the hducial region 
is dehned by: 66 GeV < mu < 116 GeV, |p^| < 2.4 and pip > 20 GeV. 

The results of our resummed calculation are shown in Fig. 0 (a) and (b). The blue-solid 
(red-dashed) histogram is the NNLL-I-NNLO (NLL-I-NLO) prediction for the qt spectrum, which 
is normalized to the cross section in the hducial region, and the points are the data with the 
corresponding experimental errors. The inset plot shows the high-g^ region. To facilitate the 
comparison between the data and the perturbative calculation we consider their ratio with respect 
to a reference theoretical result. We choose the NNLL-I-NNLO result at central values of the 
scales {pf = fJ'R = Q = as reference result. The lower panel shows the data and the scale 

dependent NNLL-I-NNLO prediction normalized to this reference theoretical prediction. The 
scale dependence band of the perturbative calculation is computed by varying pi?, pn and Q as 
previously discussed in Sect. 13.11 we vary pp, and Q in the range mzj^ < 
with the constraints 0.5 < {pf !< 2. We see that our perturbative calculation is 
consistent with the data within the uncertainties. The scale variation bands at NLL-I-NLO and 
NNLL-I-NNLO accuracy overlap. Moreover, in going from NLL-I-NLO to NNLL-I-NNLO accuracy 
the perturbative uncertainty is reduced and the agreement between experimental data and theory 
prediction is improved. The perturbative uncertainty at NNLL-I-NNLO accuracy is about ±10% 
at the peak, it decreases to about ±4% at qp ~ 10 GeV, and it increases again to about ±10% 
at Qt = 40 GeV. The comparison between our theoretical prediction and the GMS and ATLAS 
data is qualitatively similar, the main difference being that, due to the larger data sample, the 
experimental errors in the ATLAS analysis are significantly smaller. 

We add a comment on the large-gr region (see inset plots of Fig. [5|), where the cross section 
is dominated by the hxed-order contribution. For very large qp, i.e. qp 3> mz, the physical hard 
scale of the process is of the order of qp and not of the order of mz, and a sensible scale choice 
is Pi? ~ ~ qp. Therefore, it is not unexpected that our NNLL±NNLO calculations, which use 

hF ~ hR ~ slightly overshoot the GMS and ATLAS data in the last few high-gi- bins. The 

size of the QGD corrections evaluated with Pf ~ Pi? ~ would be smaller. Moreover, in the 
extreme region qp S> mz a resummation of enhanced large-gr perturbative terms is in principle 
required [871 EH] • 

In Fig.IHlwe consider the qp spectrum of bosons. We present a comparison of our resummed 
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Figure 5: Vector boson production at the LHC with lepton selection cuts. The NLL+NLO (red) 
and NNLL+NNLO (blue) normalized qt spectra for Zproduction are compared with the CMS 
data of Ref. fUT) (left panel) and the ATLAS data of Ref. l9^ (right panel). The scale variation 
bands are obtained as described in the text. The inset plot shows the ratio of the data and of the 
scale dependent NNLL+NNLO result with respect to the NNLL+NNLO result at central values of 
the scales. 


results with the pp ^ W -+ H data collected by the ATLAS Collaboration |93] with an integrated 
luminosity of 31 pb“^ at ^/s = 7 TeV. The hducial region is dehned as follows: the charged lepton 
has transverse momentum pip > 20 GeV and pseudorapidity |? 7 *| < 2.4, the missing transverse 
energy is plf > 25 GeV, and transverse mass = ■\/2p^py(l — cos(0^ — 0^)) is constrained in 
the region > 40 GeV. We recall that, because of the presence of the neutrino in the hnal state, 
the qt of the W has to be reconstructed through the transverse energy of the hadronic recoil, which 
has a poorer experimental resolution than that of the lepton momentum. In the small qt region, 
the bin sizes of the experimental data are rather large, with only four bins in the region with 
qt < 55 GeV. For this reason in Fig. [6] we focus on the large qr region 55 GeV < qt < 300 GeV, 
while the small qt region is shown in the inset plot. The lower panel of Fig. [6], which covers the 
entire qt region of the data, presents the ratio of both data and theoretical results with respect 
to the reference theoretical result. This ratio and the scale variation bands are computed exactly 
in the same manner as in the case of Fig. [5l Looking at the ratio plot in the lower panel, we 
see that our NNLL+NNLO calculation describes the W production data within the perturbative 
uncertainties. The NNLL+NNLO perturbative uncertainty is about ±8% at the peak, it decreases 
to about ±4% at ~ 15 GeV, and it increases again to about ±15% at qt = 50 GeV. 

In Sect. l3.1l and in the hrst part of this Section, we have examined vector boson distributions 
(without and with the application of acceptance cuts) and we have computed and studied the 
effects that are produced by the all-order resummation of large logarithmically-enhanced terms 
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Figure 6: Vector boson production at the LHC with lepton selection cuts. The NLL+NLO (red) 
and NNLL+NNLO (blue) normalized qt spectra for production are compared with the ATLAS 
data of Ref. f9&) . The ratio in the lower panel and the scale variation bands are obtained as in 
Fig. 0 


at small values of qt- Our related calculations are performed at complete NNLL+NNLO (and 
NLL+NLO) accuracy. In the following part of this Section, we consider other observables that 
are related to the qt distributions but in which hxed values of qt are not directly measured. 
These observables are inclusive over qt within certain qt ranges. Since the bulk of the vector 
boson cross section is produced at small values of gr, if the observable (indirectly) probes the 
detailed shape of the production cross section in the small-g^ region, the observable itself can 
be very sensitive to high-order radiative corrections and to the qt resummation effects that we 
can explicitly compute. This reasoning illustrates and justihes the physical (and quantitative) 
relevance of qt resummation for other g'r-related observables. In the second part of this Section 
we study the quantitative impact of qx resummation on some observables. 

At the formal level, our study of other observables implies that we are resumming high-order 
logarithmic corrections (in case they are present) that appear in the computation of those ob¬ 
servables. Strictly speaking, this resummation has to be performed on an observable-dependent 
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basis (see, e.g., Ref. [HS])- Therefore, our observable-independent treatment (based on transverse- 
momentum resummation) cannot guarantee that we formally achieve exact NNLL-I-NNLO ac¬ 
curacy for all these observables. Nonetheless we are able to correctly take into account all the 
leading-logarithmic contributions, all the complete (with and without logarithmic enhancement) 
perturbative terms up to the NNLO levei«], and a substantial part of subleading logarithmic terms 
beyond the NNLO accuracy. 

This statement about resummation is a consequence of the following discussion. The observable- 
dependent logarithmic terms (in case they are present) are due to multiple radiation of soft and 
collinear partons in the inclusive hnal-state: these logarithmic corrections are computed by approx¬ 
imating the QCD scattering amplitudes in the soft and collinear limits and, then, by integrating 
the hnal-state QCD radiation over the corresponding phase space with appropriate (observable- 
dependent) kinematical approximations. In our transverse-momentum resummation procedure we 
correctly take into account the NNLL dynamics (the behaviour of the QCD scattering amplitudes 
in the soft and collinear limits) of soft and collinear radiation, and we treat the phase space of the 
hnal-state QCD radiation with consistent kinematical approximations that are specihc of the qt 
spectrum. However, the observable-dependent kinematical approximations can only diher beyond 
the leading-logarithmic level (to leading-logarithmic level, a strong-ordering approximation in the 
energy/angle of the emitted partons is sufficient), and these diherences do not spoil the leading- 
logarithmic accuracy of our resummed calculation. In this respect, it is important to remark the 
role of the qt recoil (see Sect. |2]and Appendix on the kinematics of the produced (observed) 
lepton pair. We treat the qx recoil in a kinematically consistent way (though it necessarily involves 
non logarithmic approximations that are uniformly of 0 {qT/M) throughout the small-gy region), 
and such a treatment is necessary to correctly correlate the dynamical gr resummation effect with 
the ensuing dependence of the measured (computed) observable. 

In summary, the application of our q^ resummed calculations to the computation of other ob¬ 
servables is physically (and, thus, quantitatively) and formally (as we have just discussed) justihed. 
A detailed specihcation of the subleading-logarithmic accuracy of the gr resummed calculation at 
the formal (analytical) level requires (and deserves) observable-dependent investigations, which 
can be performed in future studies. 

Among other observables, we hrst consider the measuremenili] of the cf* distribution from 
pp —)■ Z/ 7 * —)■ data at ^/s = 7 TeV as reported by the ATLAS Collaboration [9l] with an 

integrated luminosity of 4.6 fb“^. The (f* observable is dehned as (f* = tan( 7 r /2 — A0/2) sin(0*), 
where A(/) is the azimuthal angle between the leptons and the angle 9* is dehned by cos 9* = 
tanh((? 7 ^^ — rf )/ 2 ) where rf"^ (rf ) is the rapidity of the positively (negatively) charged lepton. 
The cuts that dehne the hducial region are those of the ATLAS analysis of the gr spectrum: 
66 GeV < mil < 116 GeV, p^ > 20 GeV and |p*| < 2.4. 

The 0* variable at small values of 0* is correlated to g^ and, therefore, it is strongly sen¬ 
sitive to gr resummation effects. A detailed discussion on the relation between 0* and qr is 

Wpor observables that are inclusive over the region that includes qt = 0, the NNLO accuracy is achieved through 
our detailed matching procedure (see Sect. [2]) with the fixed-order calculation. 

lAn analogous measurement of the 0* distribution at the LHC was reported by the LHCb Collaboration [5^ 
with an integrated luminosity of 0.94 fb~^. In the small-0* region, the bin sizes of the LHCb measurement are 
rather large (with respect to those of the ATLAS measurement [94]), with only two (four) bins in the region 
0 * < 0.1 ( 0 * < 0 . 2 ). 
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Figure 7: The NLL+NLO (red) and NNLL+NNLO (blue) normalized 0* distribution for Z/ 7 * 
production at the LHC is compared with the ATLAS data of Ref. The NLL+NLO and 

NNLL+NNLO central results are computed at the scales fiR = pip = Q = mzj^. The ratio in the 
lower panel and the scale variation bands are obtained as in Fig. 0 


presented in Ref. |96], where the resummation of the Inc/)* terms is carried out in analytic form 
up to NNLL+(!7 (q;|) accuracjlil, and it turns out to be strictly related and very similar to qp 
resummation. Ensuing phenomenological studies are presented in Refs. [671 l 68 ]. 

In Fig. [7] we report the ATLAS data of the 0* distribution (normalized to the measured cross 
section in the hducial region) and the comparison with the results of our resummed calculation. 
The NLL+NLO and NNLL+NNLO central results are computed at the scales Pr = Pf = Q = 
mz/2. The scale variation bands at NLL+NLO (red) and NNLL+NNLO (blue) accuracy and 

t+he analytical treatment of Ref. [96] does not reach complete NNLO accuracy at small values of 0* since the 
analogue of the vector boson coefficient 'Hy'^ in Eq. ® is not included in the calculation. An approximated form 
of Ry^ is included in the calculation of Ref. [55] . 
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the reference NNLL+NNLO result for the ratio in the lower panel of Fig. [3 are computed as in 
Figs. |5] and [6l We observe that the scale variation bands at the two subsequent orders overlap, 
and that the NNLL+NNLO perturbative uncertainty is substantially smaller than the NLL+NLO 
one. The NNLL+NNLO result is consistent with the data within the uncertainties in both the 
small-0* and large-0* regions (the large-0* region is shown in the inset plot). The NNLL+NNLO 
perturbative uncertainty is about ±10% for 0* < 0.01, it decreases to about ±5% at 0* ~ 0.05, 
and it increases again to about ±10% at 0* ~ 0.2. 

We add a comment on the results that we have shown in Figs. [5H71 We recall that all the 
results presented in this Section are obtained in a purely perturbative framework. In Sect. 13.11 
we have discussed the possible impact of the inclusion of a NP form factor, and we have seen 
(Fig. [3]) that NP effects should lead to a deformation of the perturbative result that is well within 
the scale variation uncertainties of the NNLL+NNLO calculation. In Figs. [3H7]we observe that 
all the resummed perturbative predictions are consistent with the data within our estimation 
of perturbative uncertainties. Owing to the agreement between the theoretical NNLL+NNLO 
predictions and the experimental data in the very small q'r/0* region, we cannot draw any precise 
quantitative conclusion about the dehnite size of NP effects in the ^’/y*, and 0* distributions 
that we have considered. We can only conclude that NP effects have to be small in order not to 
spoil the agreement between the data and the corresponding NNLL+NNLO results in Figs. [5H71 

We conclude this Section by considering other observables. We study the impact of qt resum¬ 
mation on the kinematical distributions that are relevant for the measurement of the W mass. 
We consider pp —)■ W~ —)■ l~vi with i/s = 7 TeV and we apply the following selection cuts: the 
charged lepton has transverse momentum plp > 30 GeV and rapidity \p^\ < 2.4, the missing trans¬ 
verse momentum is > 30 GeV, and the transverse mass rriT has rriT > 60 GeV. We also apply a 
cut, p^ < 30 GeV, on the transverse momentum p^ of the W boson (lepton pair). The results of 
our calculation of the mrp distribution and of the lepton momentum distributions are presented in 
Fig. [Eland [9], respectively. The reference scale choice of the calculation is pp = = Q = mw/2. 

In both hgures we present the results of the hxed-order calculation at LO (cyan dotted), NLO 
(green solid) and NNLO (black dot-dashed) accuracy and we compare them with the results of the 
qt resummed calculation at NLL+NLO (red dashed) and NNLL+NNLO (blue solid) accuracy. 
The lower panels show the ratio between the various results and the NNLL+NNLO result (the 
ratio LO/(NNLL+NNLO) is not reported in the lower panels). 

The rriT distribution in the range mp < 90 GeV is presented in Fig. [HI We can consider two 
regions: the large-m^ region, around mp ~ mw (we recall that we use rriw = 80.385 GeV), and the 
smaXi-rrip region. In the large-m'^ region, mp > 70 GeV, we see that the perturbative prediction is 
extremely stable against radiative corrections, and the stability is present both in going from NLO 
to NNLO accuracy and with inclusion of resummation. This is a consequence of the well known 
fact that the transverse mass is weakly sensitive to the transverse momentum of the W boson. 
Formally, the mp distribution has no logarithmic corrections of the type ln(|mr — mwl/^w), and 
our qp resummed calculation does not spoil the stability of the hxed-order expansion. On the 
contrary, in the small-mr region, we observe that the hxed-order predictions become unreliable. 
The LO distribution is large at mp = 60 GeV, and both the NLO and NNLO distributions 
become negative at mp ~ 60 GeV. This (mis-)behaviour is due to the fact that the constraints 
Pp > 30 GeV and pp > 30 GeV produce an unphysical boundary (and a stepwise behaviour) of 
the mp distribution at mp = mpstep = 60 GeV in the LO calculation. The boundary is due to 
the LO kinematics Pp + Pp = Qt = and it disappears at higher orders since % 0. The LO 
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Figure 8: Effect of qt resummation on the transverse-mass (mx) distribution for pp —)■ W~ —)■ l~t)i 
production at the LHC. Comparison of results of the fixed-order calculation at LO (cyan dotted), 
NLO (green solid) and NNLO (black dot-dashed) with the resummed calculation at NLL-hNLO 
(red dashed) and NNLL-hNNLO (blue solid) accuracy. The lower panel shows the ratio between 
the various results (excluding the LO result) and the NNLL-hNNLO result. 


boundary induces (integrable) logarithmic singularities of the type ln(l — mxstep/mTY at NLO 
and beyond [98]. These logarithmic terms are resummed to all order by qx resummation, and the 
singularities are absent in the resummed prediction [98], which is well behaved at the LO boundary 
mx = rnxstep- We also note that the differences between the NLL+NLO and NNLL+NNLO results 
are small at mx ~ 60 GeV. 

In Figs, ini (a) and (b) we present the pip and pf. distributions, respectively. In the limit in 
which the W boson is produced on shell, these distributions have an LO kinematical boundary at 
mw/‘^- The hnite width of the W boson (partially) smears this effect: at LO both the pip and plf 
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(a) 



(b) 


Figure 9: Effect of qt resummation for pp —)■ W~ —)■ l~ni production at the LHC: (a) lepton pt 
distribution and (b) missing px distribution. The fixed-order and resummed results are denoted as 
in Fig. 0. 


distributions are strongly peaked at /2 (Jacobian peak) and quickly drop for px ^mvy/2. The 
almost stepwise behaviour of the LO distribution produces large radiative corrections at NLO and 
beyond (in the limit in which the W boson is produced on shell, these large corrections would 
be integrable logarithmic singularities at each perturbative order [HE]). The NLO and NNLO 
distributions indeed display an unphysical peak at px ~ 42 GeV, which is an artifact of such 
large corrections (singularities in the on-shell limit). The resummed predictions at NLLj-NLO 
and NNLL-I-NNLO accuracy are free of such instabilities and display a smooth short/der behaviour 
around the LO boundary for on-shell production. The perturbative instabilities of the hxed-order 
calculation at small values of px {px ~ 30 GeV and plf ~ 30 GeV) are analogous to those that 
we have previously discussed in the case of the mx distribution in the region mx ~ 60 GeV (see 
Fig. 0). In the case of the px distributions, it is the constraint mx > 60 GeV that produces 
the LO boundaries at Px = Pt = 30 GeV, an LO stepwise behaviour and ensuing instabilities 
at each subsequent perturbative order. The resummed calculation is perturbatively stable in the 
small-pr region, and the differences between the NLLj-NLO and NNLL-I-NNLO results are small 
throughout the entire region with px ^ 45 GeV. In the large-pr region {px ^ 45 GeV) both the pip 
and plf distributions display radiative corrections that are relatively large. This is not unexpected 
since in this region of transverse momenta the NLO calculation is essentially the hrst perturbative 
order at which both the pf and the plf distributions are non vanishing (in the on-shell limit, the 
0{as) and C>(tt|) result would be an LO and an NLO prediction, respectively). 






























4 Summary 


In this paper we have considered the transverse-momentnm (qx) distribntion of DY high-mass 
lepton pairs prodnced, via Z and W bosons decay, in hadronic collisions. We have presented 
a pertnrbative QCD stndy based on transverse-momentnm resnmmation np to NNLL accnracy. 
We have combined small-gT resnmmation with the known 0{a‘^) fixed-order resnlt at small, in¬ 
termediate and large valnes of qt- 

We have followed the resnmmation formalism developed in Refs. [25l |26l |27] to implement 
transverse-momentnm resnmmation and the matching with the result at 0{a‘g). In particular, 
our calculation includes the complete NNLO contributions at small values of qt (i.e., in any regions 
that include qt = 0) and it exactly reproduces the complete NNLO total cross section after inte¬ 
gration over Qt- This leads to theoretical predictions with a controllable and uniform perturbative 
accuracy over the region from small up to large values of qt- At large values of qt, the predictivity 
of small-gT resnmmation is superseded by that of the customary fixed-order expansion, and our 
resummed calculation can be smoothly joined onto the 0(a|) calculation. The resummed calcu¬ 
lation can be systematically expanded at various orders of logarithmic accuracy (e.g., NLL-I-NLO 
and NNLL-I-NNLO accuracy), and its theoretical uncertainties due to uncalculated higher-order 
QCD corrections can be studied by comparing the results at two subsequent orders and by per¬ 
forming systematic studies on factorization, renormalization and resnmmation scale dependence. 
We have performed such a study for the case of vector boson production at LHC energies, and 
we have briefly illustrated the uncertainties due to parton densities and the possible impact of 
non-perturbative effects. 

In the present paper we have extended the resummed calculation presented in Ref. [10] for Z/y* 
production by considering also production and by including the leptonic decay of the vector 
boson with the corresponding spin correlations, the hnite-width effects and the full dependence 
on the final-state leptonic variables. We have compared our resummed results for ^’/y* and W 
production with some of the available data of the ATLAS and CMS experiments at the LHC, 
applying the same kinematical cuts on final state leptons that are considered in the experimental 
analyses. We find that the data are well described by our predictions within the perturbative 
uncertainties. We have also considered the impact of transverse-momentnm resnmmation on 
observables, which are different from the vector boson qx, that depend on the lepton kinematical 
variables. In particular, we have studied the (p* distribution in 2’/y* production and the leptonic 
transverse-momentnm, the missing transverse-momentnm and the transverse-mass distributions 
in W production. 

Our calculation is implemented in the parton-level Monte Carlo numerical code DYRes which 
allows the user to apply arbitrary kinematical cuts on the vector boson and the final-state leptons, 
and to compute the corresponding relevant distributions in the form of bin histograms. These 
features make our program a useful tool for DY studies at the Tevatron and the LHC. A version 
of the DYRes code is publicly available. 

The production and decay mechanisms of the vector boson are dynamically correlated by the 
non vanishing spin of the vector boson. The inclusion of the lepton decay (with the spin correla¬ 
tions and the full dependence on the kinematical variables of the two leptons) in the resummed 
calculation requires a general theoretical discussion on the qx recoil due to the transverse mo- 
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mentum of the vector boson. This discussion is not limited to the specihc case of vector boson 
production. We have presented a general and explicit procedure to treat the qt recoil. The 
procedure is directly applicable to qt resummed calculations for production processes of generic 
high-mass systems in hadron collisions. 
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A Appendix: Lepton angular distribution and qx recoil in 
transverse-momentum resummation 


This Appendix is devoted to the g^-recoil issue that we have introduced and illustrated in Sect. [2] 
(see Eqs. and accompanying comments). To our knowledge the issue has not received 

much attention in the previous literature on transverse-momentum resummation. We present a 
detailed discussion of the issue and a general, explicit and consistent procedure to implement 
the Qt recoil in transverse-momentum resummation. Our procedure explicitly exhibits the degree 
of freedom involved in the implementation of the qt recoil and, moreover, it gives an explicit 
formal parametrization of the ensuing ambiguities. The procedure is straightforwardly applicable 
to implement the qx recoil (and, possibly, estimate related uncertainties) in calculations based on 
transverse-momentum resummation. 


The gT-recoil issue is not specihc of the lepton angular distribution for vector boson decay, but 
it regards transverse-momentum resummation for generic production processes. For simplicity of 
presentation, in the following we consider in detail vector boson production and the DY process. 
Then we discuss the generalization to generic processes. 


We begin our discussion by considering the computation of the DY multidifferential cross 
section in Eq. (j3]) at the LO in perturbative QGD. At this order the hadronic cross section (and the 
corresponding partonic cross section) is directly and exactly (i.e., with no small-gr approximation) 
proportional to the Born level angular distribution /dfl in Eq. ([5]). We have 
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For the purpose of our general discussion of qt recoil, we write the lepton angular distribution 
in the following form: 
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Note that, following the general notation of Eq. ([3]), we have not specihed the actual dehnition 
of the angular variables f2, and we have written the left-hand side of Eq. fl20|) in a Lorentz 
invariant form. The relation fl20|) is written in the form of a proportionality relation: the additional 
proportionality factors that are not explicitly denoted in the right-hand side are not relevant for 
our following discussion (in particular, they are independent of the lepton momenta {^ 3 ,^ 4 } and, 
thus, of n). The factor is the square of the Born level scattering amplitude 

for the partonic process 

fll(fcl) + ^2(^2) ^3(^3) + ^4(^4) ; ( 21 ) 

where ki {i = 1 , 2 ) is the momentum of the colliding parton a* from the initial-state hadron hi{Pi) 
(see Eq. ([1])), with the kinematics 


ki + k 2 = q , k^ = 0 {i = 1 , 2 ) 


( 22 ) 


In our specihc case of vector boson production, the Born level partonic process is the qq annihi¬ 
lation process qjq^ V ^ I 3 I 4 (i.e., { 01 , 02 } = I?/, ?/'})• All the other factors in the right-hand 
side of the relation fl 20 l) are related to the kinematical phase space of the hnal-state leptons and, 
in particular, they enforce the kinematical constraint q = Ps + Pi- 


The LO calculation of the cross section kinematically relates the parton and hadron momenta 
ki and P*. In particular, at the LO we have qx = 0 and, specihcally, the LO value ki of fhe 
parton momentum i^ 


^^LO) = > ^l(LO) = ^2^2" , 


(LO) 


with 


M 


Xi = 


X2 = 


M e-y 

7 ^ 


( 23 ) 


( 24 ) 


where (see Sect. |2]) M and y are the invariant mass and the rapidity of the lepton pair and ^/s is 
the hadronic centre-of-mass energy. Inserting the LO expression fl23|l of the parton momenta ki in 


Eq. 


, the LO lepton angular distribution 




LO 


in Eq. flT^ is uniquely specihed. 


Higher-order perturbative contributions produce logarithmically-enhanced (‘singular’) terms 
at small qx that can be resummed to all orders, leading to the resummation factor W in Eq. (15|) . 
These logarithmic terms are due to multiple radiation of soft and collinear partons, and this soft 
and collinear radiation is factorized [29] with respect to the Born level amplitude of 

Eq. fl20ll . As a consequence, after qx resummation the angular distribution of the decaying leptons 
is still given by the Born level function /dO. in Eq. fl20l) . and this function thus appears as 
a multiplicative factor in front of the resummation factor W of the resummed component of the 
vector boson qx cross section (see Eq. ([5|)). Strictly speaking |29|, in the limit qx ^ M that is 
relevant for resummation, the angular distribution can be expressed in terms of the LO distribution 
[da^^ydn]^^ in Eq. (IT^ . namely the expression (j2nil with the LO kinematics of Eqs. (j23l) and 
(12T|) . which in particular has qx = 0. Indeed, after soft/collinear factorization and resummation, 
any residual dynamical effect on the process in Eq. (!?!]) (and on and da^^'* /dfi) is due to 
hard-parton radiation. Hard radiation produces 0{qx/M) corrections that lead to non-singular 
contributions if qx -C M: these corrections can be formally approximated by their limiting 

§The kinematical variables xt in Eqs. (ESD and (1^ and the kinematical variable zi in Eq. (1^ should not be 
confused with the integration variables Xi and Zi used in Sect. [2] (we use the same symbols for both set of variables). 
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behaviour as g-r —^ 0 and, thus, neglected in the computation of the resummed component (see 
Eqs. dH) and ffTT|) i and included in the hnite component (see Eq. (jl])). 

Neglecting these 0{qT/M) corrections is a perfectly suitable procedure for the resummed cal¬ 
culation of the vector boson qt cross section (see Eq. (inj). However, performing the resummation 
at hxed lepton momenta, the momentum of the vector boson must be fully specihed by the lepton 
momenta and, in particular, qx = pxz + Pt^ is not vanishing. The resummation factor W (see 
Eq. ([5])) produces a smearing of the LO distribution of Eq. flT^ and hnite values of qx'- 

to avoid unphysical results (e.g., events with qx ^ 0 and px^ + Pt^ = 0) the factor /dfl in 
Eq. ([5]) cannot be the LO angular distribution [dd ^^^(which has Pxs + Pt^ = 0) in Eq. ([19]). 
In other words, the non-vanishing value of qx has to be distributed between the two lepton mo¬ 
menta and this leads to the gT-recoil issue that we have illustrated in Sect. [2] (see Eqs. f[T2D - f[TTl) 
and accompanying comments). The resummed calculation requires the specihcation of a gy-recoil 
prescription that has to be consistent (and physically sensible), although this can be done in many 
(inhnitely many) different ways. 

Actual resummed calculations performed in the literature do not mention the g^-recoil issue. 
The calculations of Refs. [SSIEHIET] directly refer to the use of the Collins-Soper (CS) rest frame 
ra- The procedure to compute the factor da^^^dfl in Eq. (ED is as follows. The lepton angular 
variables Li are specihed to be the polar and azimuthal angles {O'cs^ 4>cs} leptons in 

the CS rest frame. The LO distribution in Eq. (1191) is then expressed in terms of 

{0'x:g, (j)cs} (since the LO distribution has qx = 0, in this case {O'^g, (j)cs} exactly coincide with 
the lepton scattering angles in the centre-of-mass frame of the LO colliding parton momenta 
in Eq. (1221) i and this leads to an unambiguously dehned angular function 
(this function is actually independent of (j^^g) that is used to dehne (see Eq. ([T2D i the angular 
distribution dd^^ydfl of the resummed component of the cross section (see Eq. ([5D). This is a 
perfectly dehned and consistent procedure, but it hides the actual implementation of 0{qx/M) 
corrections through an implicit prescription for the qx recoil: the dehnition of the CS rest frame 
is qx dependent and a qx dependence is introduced by identifying/equating the angles {O^g, 0cs} 
of the LO and resummed calculations (additional comment on this are presented in a paragraph 
after Eq. (132D L 

Here we explicitly present a consistent gr-recoil procedure and an entire class of gr-recoil 
prescriptions. Our viewpoint is as follows: the non-vanishing value of qx of dynamical origin that 
is produced by resummation leads to a gy-recoil that can be ‘kinematically absorbed’ □ by the 
momenta fci and ^2 of the colliding partons of the underlying hard-scattering process (see Eq. fl2TD L 
As specihed below, there are inhnitely-many ways of implementing this kinematical recoil on the 
colliding partons in a consistent manner (i.e., without modifying the logarithmically-enhanced 
perturbative terms at small qx)'- they diher by corrections that are of 0{qx/M) order-by-order 
in the perturbative expansion (after having matched the resummed calculation with the complete 
N^LO calculation, as in Eqs. (jlD and (ITUD . these corrections start to contribute at the N^+^LO 
level). 

^The qx recoil issue does not arise in the context of transverse-momentum [kx) factorization |80j for high-energy 
(small-a;) hard-scattering processes. In this formulation the qx recoil is dynamically (and uniquely) embedded in 
the factorization formula. The parton densities of the colliding hadrons are kx dependent and the hard-scattering 
colliding partons have ensuing non-vanishing transverse momenta kix that enter as integration variables in the 
factorization formula. 
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According to our procedure, the lepton angular distribution to be used in the re- 

summed calculation (see Eq. ([5])) is exactly given by the expression in Eq. fl20|) . The phase space 
factor in the right-hand side of Eq. fl20|) is directly given in terms of the physical (measured) lepton 
momenta ps and p 4 (with ps + Pa = <?)• The momentum ki (then, k 2 = q — ki) to be used to 
compute the Born level scattering amplitude in Eq. (l20ll is given by the following parametrization: 


K = zi 




2q-Pi 


Pi' + k^r, + 


^ Q- Pi pM 

M^Pa-P 2 ' 


{kirpki'j'^ — ^ir) ’ 


(25) 


where 


~ , (Af| ^ AP + #) , (26) 

and k^rp is a two-dimensional vector that is transverse to both P^ and (i.e., kix lies in the Qt 
plane) and that fulhls the following constraints: 


kij' —)■ 0 if (jx —t 0 , 
M'^ ‘2,Qx ■ kix ^ 2 Mx l^irl 


(27) 

(28) 


We note that, following the dehnition in Eqs. fl25|) and fl26|) . /cf and k^ are well dehned ‘physical’ 
parton momenta: they fulhl the kinematics in Eq. fl22l) and they have positive dehnite energies, 
> 0 and fc® > 0 (the constraint in Eq. fl25D guarantees that the four-momentum k'i has positive 
dehnite energy and, then, k^ > ^ follows from q^ > 0). Therefore the scattering amplitude 
^ 2 ;P 3 ,P 4 ) in Eq. (1201) is well dehned and unambiguously computable. Moreover, due 
to Eq. fl271) . the parton momentum ki in Eq. fl25ll coincides with its LO expression fl2^ if qx = 0. 
We also note that fcf is invariant under longitudinal boosts of the hadronic centre-of-mass frame, 
provided k^^p is boost invariant. 


At hxed values of q^, Pf and P 2 , Eqs. (I2^ -(I28 |) give the most general expression of ki that 
respects the Born level kinematics in Eqs. and fl2^ and the LO kinematics in Eqs. fl2^ and 
(EH). This expression is parametrized by the arbitrary (though constrained) transverse-momentum 
vector kix- By choosing diherent values of kix, we can obtain an entire class of consistent g^-recoil 
prescriptions. For example, two ‘obvious’ possible choices are as follows: 


A) set kiT = gr/2 (and thus = gr/2): 


from Eq. fl26|) we obtain 


and we have 


Mx + M q ■ Pi q ■ P 2 
Mr M2 Pi ■ P 2 


ei 




Mr + M g ■ P 2 , 1 , Mr - M 

2Mr Pi ■ P 2 2Mr 

Mr — M q ■ P 2 p /2 , ^ /i , Mr -|- M 
2Mr Pi ■ P 2 ^ ^ 2 2Mr 


q • Pi 
Pi-P2 
q • Pi 
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B) set kix = 0 (and thus k 2 x = Qx)- 
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from Eq. fl26p we obtain Zi = 1 and we have 


k^i 


2q-Pi 


TDt^ 

1 ) 


2q-Pi 


P^ . 


We also note that, after integration over the lepton angnlar variables Cl, we consistently obtain 
the Born level total cross section of the resnmmation formula flTTl) . Indeed, after 

the Cl integration of Eq. fl20l) . the result does no longer depend on the lepton momenta and, 
since it is a Lorentz invariant quantity, the result can only depends on the invariant (fci + /c 2 )^ = 
2ki ■ k 2 = q^ = M^, which is independent of ki^. In other words, the dependence on the arbitrary 
parameter kiT completely cancels in lepton-inclusive observables. 

Using our g^-recoil procedure, we can compute the corresponding lepton angular function 
of Eq. ([12]). This function is the product of two factors. One factor is a purely kine- 
matical origin (it derives from the phase space factor in the right-hand side of Eq. fl20|) ). and it 
depends on the specihcation of the angular variables f2. The other factor, denoted as in 
the following (for simplicity we omit the subscript g/g/' hU), has a dynamical origin and it 
depends on the Born level factor \M^^\ki, k 2 ',P 3 ,P 4 )\‘^ in the right-hand side of Eq. fl20|) . Since 
\M^^\ki, k 2 ;p 3 ,Pi)\‘^ is a Lorentz invariant scalar quantity and the momenta {ki, k 2 ;p 3 ,P 4 } are 
constrained by momentum conservation, can only depend on the dimensionless variable 

4^1 • P3/{2ki ■ /C 2 ) = 4A;i • Considering the centre-of-mass frame of ki and k 2 , we have 

Aki ■ ps/M^ = 1 — cos 0 ( 3 , where is the scattering angle between ki and P 3 . In other words, 

p{D) ^ 

and 6^(3 is the lepton scattering angle in a particular rest frame of the vector 
boson momentum q^ (the centre-of-mass frame of ki and k 2 ). Our gr-recoil procedure can thus be 
reinterpreted in terms of generation of lepton-pair events. Considering a dehnite (with respect to 
the hadronic collision frame) rest frame of the vector boson momentum, the lepton-pair event and 
the individual lepton momenta are generated in that frame according to the corresponding Born 
level angular distribution; then the lepton pair distribution is boosted to the hadronic collision 
frame through the corresponding Lorentz transformation. Since there is an inhnite numbers of 
vector boson rest frames, this event-generation procedure has an inhnite degree of arbitrariness. 
Applying a three-dimensional rotation to a vector boson rest frame, we obtain another vector 
boson rest frame and, thus, the inhnite numbers of vector boson rest frames depends on the two 
scalar parameters of the three-dimensional rotation. Accordingly, our g^-recoil procedure depends 
on the arbitrary two-dimensional vector kix, namely on two parameters (the magnitude l^irl and 
the azimuthal angle of kir). In other words, the relation between the LO momenta fcf(LO) 

Eq. (12^ and the g^-recoiled momenta ki obtained through Eqs. fl2SD - (l25D can be reinterpreted as 
a Lorentz transformation of the colliding parton momenta from the hadronic collision frame to a 
specihed vector boson rest frame. This interpretation directly relates our g^-recoil procedure with 
the specihc CS frame procedure (as already mentioned and described in the initial part of this 
Appendix) that is directly used in other resummed calculations [55] EG] |57| . It can be explicitly 
checked that the CS frame procedure used to specify the lepton angular distribution /dCl in 
the resummed calculation of Refs. [SEl ESI EZ] corresponds to the choice kix = k^r = Qt/2 (see 
Eqs. (130|) and ([3T])) within our class of gT-recoil prescriptions. 

Owing to our explicit parametrization and implementation of the gr-recoil procedure, the 
quantitative effects of various gT-recoil prescriptions can be directly investigated in applications 
of the numerical program DYRes (and of other resummed calculations). Obviously (as discussed in 
Sect. ED, different g^-recoil prescriptions have no effects on quantities that are fully inclusive over 
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the leptonic variables f2. In general, our expectations are as follows. We expect that the quan¬ 
titative differences produced by various qT-recoil prescriptions are small for lepton non-inclusive 
observables that are mostly sensitive to either the small-g^ region (in this region the gy-recoil 
effects are non-singular and thus subdominant with respect to the singular logarithmic contri¬ 
butions) or the liigh-gT region (in this region the gr-recoil effects are suppressed by the smooth 
switching procedure of Eqs. (1T5|1 - (|T7I) ), while relatively larger differences can appear in case of 
sensitivity to the region of intermediate values of qt- Moreover, these quantitative differences are 
expected to decrease in going from NLL-I-NLO to NNLL-I-NNLO accuracy, since the non-vanishing 
^T-recoil effects start to formally contribute at the level in the N^LL-|-N^LO calculation. 

In Sect. 13.21 we have presented our quantitative results obtained with the DYRes code. As stated 
at the beginning of Sect. [3l these results are obtained (analogously to those in Refs. [55l EHl [57]) 
by computing the Born level angular distribution in the CS rest frame, i.e. by setting 

kiT = kax = Qt/‘^ in the actual implementation of our g^-recoil procedure (this corresponds to 
use the prescription A in Eqs. fl29|) - fl^ h Setting kix ^ qT/‘^, we have also considered other 
variants of the qx-Tecoil prescriptions and we have examined the quantitative differences that are 
produced on the observables that are examined in Sect. 13.21 (i.e., the observables in Figs. [5HH]). We 
have found that various g^-recoil prescriptions produce differences that are in agreement with our 
general expectations and, in particular, at NNLL-I-NNLO accuracy these differences lead to small 
quantitative effects: typically, the effects are much smaller than the scale-variation uncertainties 
(estimated as in Sect. 13.2p . For instance, comparing the qr-Tecoil prescriptions A (see Eqs. fl2^ - 
fl3T]l ) and B (see Eq. fl3^ ). we obtain quantitative differences that are at most at the percent level 
(e.g., in the case of the (j)* distribution of Fig. [7] in the region 0.3 < 0* < 1, and in the case of the 
lepton-pr and missing-p^ distributions of Fig. [9]in the region 45 GeV ^pr ^ 50 GeV): these differ¬ 
ences are dehnitely smaller than the scale uncertainty (at both the NLL-I-NLO and NNLL-I-NNLO 
levels) and, especially, they are of the same size as (and, hence, hardly distinguishable from) the 
pure numerical errors of the DYRes calculation in the NNLL-I-NNLO mode. 

The ^T-recoil issue that we have discussed in this Appendix is not specihc of vector boson 
production and the ensuing leptonic decay. The issue affects qr resummed calculations for any 
process of the type hi + h 2 ^ J^{.P 3 ,P 4 ,P 5 ,...) + W (we use the same notation as in Eq. (fT|)) 
where the hnal-state high-mass system has total transverse momentum qx and the momenta 
P3,P4,P5,... of its ‘decay products’ are directly measured. Owing to the universality (process- 
independent) structure of transverse-momentum resummation [21], the g^-recoil procedure that 
we have introduced in this Appendix is directly applicable to all these processes. The only key 
difference with respect to vector boson production is that the Born level scattering amplitude 
M^^\ki,k2 ;P3,P4) in Eq. (|20D is replaced by a properly computable (all-loop) hard-virtual am¬ 
plitude M(fci, fc 2 ;p 3 ,P 4 ,P 5 ,...; as'(M^)) (see Sect. 4 in Ref. [21]), which embodies QGD virtual 
radiative corrections (M is computable as power series in as{M^)). Strictly speaking [H], the qx 
resummed cross section at small values of qx is proportional to the angular dependent distribution 
of the momenta {P 3 ,P 4 ,P 5 ,... } as computed from M a.t qx = 0 (i.e., with the LO momenta fcf(Lo) 
of Eq. fl23|l L The ensuing g^-recoil issue can be directly solved by our gr-recoil procedure. Indeed, 
the hard-virtual amplitude M{ki, k 2 ;p 3 ,P 4 ,,P 5 ,...; asiM"^)) has the same kinematical properties 
as its Born level counterpart : therefore, the qx recoil can be directly implemented by 

simply evaluating M{ki, A: 2 ;P 3 ,P 4 ,P 5 ,...; as{M^)) with the gy-recoiled momenta ki of Eqs. (I2S]) - 

([2S|). 


We add some hnal comments on spin correlations and on the specihc process of SM Higgs boson 
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production and its decay in colourless particles (e.g., H —)■ 77, H —)■ WW iulu, H ^ ZZ ^ 4£) 
|43] . The qt resummed Higgs boson cross section at fixed momenta of the decay products is 
proportional to the angular distribution as obtained (analogously to Eq. fl20|) ) by the correspond¬ 
ing Born level scattering amplitude ^ 2 ;P 3 ,P 4 ,P 5 , • • •) for the gluon fusion pro¬ 

cess gg ^ H ^ hUh ■ ■ ■ ■ Owing to the spin-0 nature of the SM Higgs boson, factorizes 
in two independent factors, ^ 2 ; <?) and {(l]PiiPAiPbi ■ ■ ■)■, for the produc¬ 

tion {gig 2 —!• H) and decay {H —)■ I 3 UI 5 ...) subprocesses of the Higgs boson. We thus have 

(g;P3,P4,P5,...)|2. Note that only depends 

on observable momenta, while k 2 ', q)\‘^ only depends on {ki -|- ^ 2 )^ = because of 

Lorentz invariance. As a consequence, our g^-recoil procedure (and its dependence on the defini¬ 
tion of ki and ^ 2 ) has no effect on the angular distribution of the Higgs boson decay products. The 
angular distribution of the resummed calculation can be directly obtained |13] by supplementing 
the Born level total cross section with the (kinematical and dynamical) Higgs boson 

decay factor. This specific example also clearly illustrates that the gT-recoil issue that we have 
introduced in Sect. |2] and discussed in this Appendix is directly related and due to the vector 
boson spin and the spin correlations between the production and decay subprocesses of the vector 
boson. This kind of relation between recoil and spin correlations is valid for generic production 
processes. 
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